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(AsQa),  arsenic  anti  site-plus-arsenic  vacancy  (AsGaVAs)  complex  as  well  as 
oxygen  complex  have  been  considered  by  many  researchers  as  the  likely  candi¬ 
dates  for  GaAs  grown  under  As-rich  or  high  arsenic  pressure  condition.  How¬ 
ever,  the  subject  is  still  highly  controversial,  and  quantitative  modelling  of 
the  EL2  electron  trap  is  needed  in  order  to  control  this  native  defect  in  both 
the  bulk  and  epitaxial  grown  GaAs. 

*'/*The  objectives  of  this  research  are  to  conduct:  (1)  A  detailed  analysis 
of  the  grown-in  defects  and  radiation  induced  defects  in  GaAs  and  other  III-V 
materials  grown  by  the  LEC,  VPE,  LPE,  and  MOCVD  techniques  under  different 
growth  and  annealing  conditions,  (2)  Theoretical  modelling  of  the  native  defect; 
for  identifying  the  physical  origins  of  the  deep-level  traps  in  GaAs  and  other 
III-V  materials,  (3)  Theoretical  and  experimental  study  of  the  potential  well 
of  electron  traps  from  analyzing  the  electric  field  enhanced  emission  rates  de¬ 
duced  from  the  nonexponential  DLTS  data,  and  (4)  Study  of  one-MeV  electron  ir¬ 
radiation  induced  deep-level  defects  in  GaAs,  AlGaAs,  and  InP  materials.  The 
main  research  accomplishments  are  summari zed^ aa  f ol towsr 

(1)  Theoretical  and  experimental  studies  of  native,  point  defects  in  GaAs  grown 
by  the  LEC,  VPE,  LEP,  and  MOCVD  techniques  under  different  growth  and 
annealing  conditions  have  been  made  in  this  research.  The  main  findings 
include:  (a)  high  purity  GaAs  material  can  be  grown  for  the  low  arsenic 
pressure  case  under  optimum  cooling  condition,  (b)  GaAs  grown  under  high 
arsenic  pressure  condition  will  produce  more  native  point  defects  than 
under  lower  arsenic  pressure  condition,  (c)  arsenic  antisite  (As Qa)  defect 
can  only  be  produced  in  VPE  and  MOCVD  GaAs  grown  under  As-rich  or  high 
arsenic  pressure  condition;  this  defect  can  not  be  produced  in  LPE  GaAs 
grown  under  low  arsenic  pressure  condition. 

(2)  A  new  defect  model  supported  by  the  experimental  data  has  been  developed 
in  this  work  to  account  for  the  physical  origins  of  EL2  electron  trap  in 
GaAs.  It  Is  shown  that  EL2  electron  trap  may  be  attributed  to  two  dif¬ 
ferent  types  of  native  defects:  One  is  identified  as  the  EL2a  (i.e., 
Er-0.83eV)  electron  trap,  and  the  other  is  designated  as  the  EL2b 
(E<j-0.76eV)  electron  trap.  The  physical  origin  for  the  EL2a  level  is  at¬ 
tributed  to  the  double-charge  arsenic  antisite  (i.e.,  AsGa++)  defect, 
whereas,  the  physical  origin  for  the  EL2b  electron  trap  is  due  to  the 
arsenic-anti  site-plus-arsenic  vacancy  complex  (i.e.,  AsGavAs)*  Based  on 
this  model,  relationship  between  the  density  of  EL2a  and  EL2b  trap  levels 
and  the  [As]/[Ga]  ratio  in  the  MOCVD  and  VPE  grown  GaAs  were  established. 
The  result  shows  the  density  of  EL2a  trap  level  is  proportional  to  the 
mole  fraction  ratio  of  (r-l)1*,  while  the  density  of  EL2  trap  level  is 
proportional  to  the  mole  fraction  ration  of  (r-1)**,  where  r  =  [As]/[Ga]  . 

(3)  A  theoretical  model  for  the  nonexponential  DLTS  response  due  to  field 
dependent  emission  rate  of  trapped  charge  has  been  developed.  A  compari¬ 
son  of  the  theoretical  calculation  of  the  nonexponential  DLTS  response 
with  the  DLTS  data  for  each  trap  level  enable  us  to  determine  the  poten¬ 
tial  well  of  the  trap  involved.  This  method  has  been  applied  to  identify 
the  physical  origins  of  the  E12a  electron  trap  in  GaAs. 

(4)  Characterization  of  the  low  energy  proton  (50  to  290  keV)  and  one-MeV 
electron  irradiation  induced  defects  in  LPE  grown  GaAs  and  LEC  grown  p- 
InP  materials  has  been  carried  out  in  this  work.  A  comparison  of  the 
deep-level  defects  produced  by  the  one-MeV  electron  irradiation  and  the 
low-energy  proton  irradiation  reveals  that  there  is  a  significant  dif¬ 
ference  In  defect  profile  as  well  as  types  of  defects  produced  in  the  GaAs 

_ LPE  layers.  The  defect  profiles  are  much  more  complicated  for  the  proton 


irradiated  GaAs  than  the  one-MeV  electron  irradiated  GaAs.  The  radia¬ 
tion  induced  deep-level  defects  are  believed  to  be  due  to  Frenkel  type 
defects  and  their  complexes.  Our  study  shows  that  low  temperature 
thermal  annealing  (i.e.,  T<300°  C)  or  periodic  annealing  process  is 
effective  in  reducing  or  eliminating  these  radiation  induced  defects  in 
both  GaAs  and  InP  materials.  The  shallow  traps  are  much  easier  to 
anneal  out  than  the  deep-level  traps. 
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I.  INTRODUCTION 


Studies  of  deep-level  defects  in  III-V  compound  semiconductors  have 
been  reported  extensively  in  the  li terature^-^] .  Speculations  on  the 
physical  origins  of  defect  levels  and  defect  complexes  in  these  materials 
are  still  very  tentative.  It  is  not  even  possible  to  say  with  any 
assurance  that  simple  gallium  vacancy  (V^)  or  arsenic  vacancy  (V^)  can 
be  related  to  a  particular  energy  level. 

In  one  respect  understanding  has  been  improved.  Gallium  arsenide 
(GaAs)  specimens  grown  by  various  techniques  such  as  liquid  phase 
encapsulation  (LEC) ,  Bridgmann,  vapor  phase  epitaxy  (VPE),  metalorganic 
chemical  vapor  deposition  (MOCVD) ,  liquid  phase  epitaxy  (LPE),  and 
molecular  beam  epitaxy  (MBE) ,  are  now  recognized  as  likely  to  have 
different  properties  in  terms  of  energy  levels  within  the  bandgap.  In  each 
of  these  techniques,  the  growth  temperature,  growth  pressure,  growth  phase, 
and  cooling  rate  are  usually  different.  For  examples,  GaAs  grown  by  the 
LPE  technique  from  a  gallium  melt  is  expected  to  be  low  in  gallium  vacancy 
defects  and  possibly  high  in  arsenic  vacancy.  The  LPE  GaAs  usually 
contains  hole  traps  with  energy  level  of  Ev+0.71eV  (B  center),  while  VPE 
GaAs  grown  in  the  As- rich  conditions  always  contains  an  electron  trap  with 
energy  of  Ec-0.83  eV  (EL2  level).  The  reason  for  observing  these  trap 
levels  in  the  LPE  or  VPE  GaAs  is  still  not  clear.  Thus,  studies  of  grown- 
in  defects  in  III-V  compound  semiconductors  are  important  tasks  to 
undertake  in  order  to  obtain  a  better  understanding  of  the  physical 
properties  of  deep-level  defects  in  III-V  materials. 

The  main  objectives  of  this  research  program  are  to  conduct:  (1)  A 
detailed  analysis  of  the  grown- in  defects  and  radiation- induced  defects  in 
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GaAs  and  other  III-V  compounds  grown  by  the  LEC,  VPE,  LPE,  and  MOCVD 
techniques  under  different  growth  and  annealing  conditions,  (2) 
theoretical  modelling  to  identify  the  physical  origins  of  the  deep-level 
traps,  in  particular  the  EL2  electron  trap  in  GaAs,  (3)  theoretical  and 
experimental  studies  to  determine  the  potential  well  of  electron  traps  from 
analyzing  the  field  enhanced  emission  rates  deduced  from  the  nonexponential 
DLTS  response  data,  and  (4)  experimental  studies  of  the  grown-in  defects 
and  one-MeV  electron  irradiation-induced  defects  in  GaAs  and  InP  under 
different  radiation  fluences  and  annealing  conditions. 

A  detailed  theoretical  and  experimental  study  of  the  grown-in  defects 
and  radiation  induced  defects  in  GaAs  and  InP  has  been  carried  out  in  this 
work  and  the  results  are  summarized  as  follows: 

(1)  Theoretical  modelling  of  native  defects  in  GaAs  has  been  developed. 
GaAs  specimen  grown  by  the  LEC,  VPE,  LPE,  and  MOCVD  techniques  under 
different  growth  and  annealing  conditions  have  been  studied  in  this 
work  in  order  to  verify  this  model.  GaAs  specimen  used  in  this 
study  include:  (i)  VPE  GaAs  epi  layers  grown  on  <100>,  <211A>,  and 
<211B>  oriented  semi -insulating  (S.I.)  Cr-doped  GaAs  substrates  with 
[Ga]/[As]  ratios  varying  from  2/1  to  6/1,  (ii)  LEC  grown  n-GaAs 
samples  which  were  used  for  study  of  the  effect  of  hydrogen  heat 
treatment  on  deep- level  defects,  (iii)  LPE  GaAs  samples  grown  under 
two  different  temperatures  (ie.,  700  and  800°C)  and  two  cooling  rates 
(0.4  and  l°C/min.),  and  (iv)  MOCVD  GaAs  epi  layers  grown  on  S.I.  GaAs 
substrate  and  S.I.  Ge  substrates.  Theoretical  modelling  of  native 
point  defects  is  described  in  chapter  V.  Experimental  results  for 
some  of  the  native  point  defects  observed  in  the  samples  cited  above 
are  described  in  chapter  VII. 


(2)  Deep-donor  trap,  commonly  known  as  EL2  center  ,  with  activation  energy 
ranging  from  Ec-0.76  to  Ec-0.83eV,  has  been  observed  in  GaAs  grown  by 
LEC,  VPE,  and  MOCVD  techniques.  The  EL2  level  acts  as  a  recombination 
center  for  lifetime  reduction.  The  physical  origin  of  this  electron 
trap  is  a  subject  of  greater  interest  in  recent  years.  Although  a 
large  number  of  papers  has  been  devoted  to  finding  the  physical 
origins  of  the  EL 2  electron  trap,  unfortunately,  none  of  these 
published  results  offered  a  convincing  and  unambiguous  explanation 
for  the  observed  EL2  trap  in  GaAs.  Our  theoretical  modelling  and 
experimental  evidence  lead  us  to  believe  that  the  EL2  center  in  GaAs 
may  be  attributed  to  two  different  types  of  native  defects:  one  is 
identified  as  the  EL2-a  electron  trap  which  has  an  activation  energy 
of  Ec-0.83  eV  and  is  due  to  the  AsGa  anti  site  defect;  the  other  is 
designated  as  the  EL2-b  electron  trap  which  has  an  activation  energy 
of  Ec-  0.76  eV  and  is  attributed  to  theAsGa-VAsdefect  complex.  A 
detailed  theoretical  and  experimental  study  of  the  EL2  electron  trap 
has  been  carried  out  in  this  work,  and  the  results  are  discussed  in 
chapter  VI. 

(3)  DLTS  signals  are  often  analyzed  assuming  that  the  capacitance 
transient  is  exponential. However,  most  of  the  DLTS  signals  are 
nonexponential  transient. [15,16]  In  ^is  Work,  we  present  the 
theoretical  analysis  of  nonexponential  capacitance  transients  due  to 
electric  field  dependent  emission  rate  of  trapped  charge.  Emission 
rate  of  trapped  charge  carriers  is  enhanced  by  the  Poole-Frenkel  and 
phonon-assisted  tunneling  effects  in  the  presence  of  an  applied 
electric  field.  Since  the  electric  field  varies  with  position  within 
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the. depletion  region  of  a  reverse  biased  p-n  junction,  the  emission 
rate  is  not  constant  within  the  same  region.  The  DLTS  response  due  to 
nonuni  formed  emission  rates  can  in  general  be  expressed  as  Sf?)= 
5 exp(-enj tx)  -£exp(-enjt2).  Based  on  the  nonexponential  capacitance 
transient,  theoretical  calculation^  of  DLTS  response  for  deep- level 
traps  in  GaAs  were  made  using  five  different  potential  wells;  namely, 
the  Coulombic  well  which  has  a  positive  charge  state  for  the  empty 
state;  the  Dirac  well,  square  well,  polarization  well,  and  dipole  well 
which  all  have  neutral  charge  state  from  different  physical  origins. 
A  comparison  of  the  theoretical  calculation  of  the  nonexponential  DLTS 
response  with  the  DLTS  data  for  each  trap  would  allow  us  to  determine 
the  potential  well  of  the  trap  involved.  This  model  was  applied  to 
the  calculations  of  potential  well  for  the  EL2a  electron  trap  and 
other  electrontraps.  The  results  will  be  discussed  in  chapter  IV. 
(4)  Characterization  of  low  energy  proton  and  one-MeV  electron  irradiation 
induced  defects  in  the  LPE  grown  GaAs  has  been  carried  out,  and  the 
results  were  described  in  details  in  the  1983  AFOSR  Annual  Technical 
Report. 

Experimental  tools  employed  in  this  study  include  the  current-voltage 
(I-V) ,  capacitance- voltage  (C-V),  thermally  stimulated  capacitance  (TSCAP) 
and  the  deep-level  transient  spectroscopy  (DLTS)  measurements.  From  these 
measurements,  one  can  determine  the  energy  levels,  density  of  defects, 
defect  profile,  and  the  capture  cross  section  for  each  trap  level. 

Chapter  II  reviews  the  possible  native  defects  and  impurity  complexes 
in  GaAs.  Chapter  III  presents  the  experimental  details.  Determination  of 


potential  well  of  deep  level  traps  using  field  enhanced  emission  rate 
analysis  of  nonexponential  DLTS  in  GaAs  will  be  described  in  chapter  IV. 
Chapter  V  depicts  the  modelling  of  the  grown-in  native  point  defects  in 
GaAs.  The  physical  origin  of  the  EL2  electron  trap  in  GaAs  is  explained  in 
chapter  VI.  Chapter  VII  described  the  study  of  grown-in  defects  vs  growth 
parameters  in  VPE,  LEC,  VPE,  and  MOCVD  grown  GaAs.  Summary  and  conclusions 
are  given  in  chapter  VIII. 


II.  REVIEW  Cf  NATIVE  BEEBCIS  AND  DBBCT  COMPLEXES  IN  GALLIOM  ARSHJIDE. 

The  number  of  possible  native  defects  in  GaAs  is  large  as  may  be  seen 
in  table  2.1.  None  of  these  native  defects  has  been  identified  with  any 
confidence.  This  is  due  to  the  fact  that  experiments  for  studying  such 
defects  tend  to  be  too  uncontrollable.  In  addition  to  the  defects  shown  in 
table  2.1r  other  types  of  defects  such  as  impurity,  complexes  of  impurity 
may  also  be  expected  in  GaAs  as  well  as  other  III-V  materials.  In  this 
chapter,  we  will  focus  our  attention  on  a  few  native  defects  listed  in 
table  2.1,  which  are  believed  to  be  related  to  the  electron  or  hole  traps 
observed  in  this  study. 

Defects  can  be  represented  by  their  chan i cal  symbols.  For  examples, 
vacancy  is  represented  by  V  and  interstitial  by  i.  Subscript  indicates  the 
lattice  site.  Thus,  vGa  denotes  the  gallium  vacancy,  Ga^  represents  the 
gallium  interstitial  site  defect.  In  addition  to  these  simple  point 
defects,  defect  complexes  (for  example:  VGa“  ASgg^VQg”)  might  be  expected 
to  form  as  the  crystal  cooling  down  from  high  growth  temperature  to  room 
temperature.  Another  type  of  defect,  namely,  the  anti  site  defect  must  also 
be  considered.  The  anti  site  defects  such  as  AsGa++,  GaAs  ,  and 
ASQg^Ga^  ,  are  believed  to  be  important  native  defects  in  GaAs. 

A  survey  of  the  literature  on  the  subject  of  defects  in  GaAs  grown  by 
the  Liquid  Encapsulation  Czockraski  (LEC),  Horizontal  Bridgman-Stockhargen 
(HB),  Vapor  Phase  Epitaxy  (VPE),  Liquid  Phase  Epitaxy  (LPE),  Molecular  Beam 
Epitaxy  (MBE) ,  and  Organometall ic  Chemical  Vapor  Deposition  (MOCVD) 
techniques  showed  that  only  a  few  electron  and  hole  traps  are  common  point 
defects  observed  in  both  the  bulk  and  epitaxial  grown  GaAs  material.  It 
should  be  noted  that  defects  observed  in  one  GaAs  specimen  may  not  be 


observed  in  other  GaAs  specimen  if  the  growth  processes  were  different.  A 
list  of  the  electron  and  hole  traps  observed  in  this  study  compared  with 
others  123-33]  are  given  respectively  in  table  2.2  and  2.3. 

The  EL2  electron  trap  (  same  as  EF10,  ET1,  EB2,  ESI  cited  in  the 
literature)  is  the  dominant  electron  trap  in  the  VPE,  MOCVD,  and  bulk  grown 
GaAs,  and  is  absent  in  the  LPE  and  MBE  grown  GaAs  epitaxial  layers.  Trap 
densities  of  EL3,  EL5,  and  EL12  levels  are  varied  by  different  growth 
conditions  ,  and  density  of  these  traps  can  be  easily  reduced  by  annealing. 
This  suggests  that  they  are  due  to  native  point  defects.  On  the  other 
hand,  EL11  is  not  affected  by  heat  treatment,  and  therefore  could  be  due  to 
an  impurity  complex.  Crystals  grown  by  different  growth  techniques  are 
expected  to  produce  different  electron  traps,  and  only  a  few  of  than  are 
common  and  independent  of  methods  of  crystal  growth.  For  example,  both  EL6 
(bulk)  and  EL7  (MBE)  levels  could  be  the  same  defect  observed  in  the  bulk 
and  MBE  grown  GaAs,  yet  it  was  not  observed  in  VPE  grown  GaAs.  EL10  (MBE) 
and  ELll  (VPE)  could  be  the  same  level,  yet  it  was  not  observed  in  bulk 
GaAs  material.  Most  of  these  deep  level  defects  are  believed  to  be  either 
due  to  the  vacancy  related  defects,  anti  site  defects  or  vacancy- impurity 
complex  defects.  For  the  hole  trap  levels  in  GaAs,  HL2,  HL5,  HL7,  HL9,  and 
HLll  may  be  related  to  native  defects,  while  the  other  trap  levels  are 
related  to  the  impurity  defects. 


Table  2.2.  Electron  traps  in  n-GaAs 


Level 

Activation 
energy  (eV) 

Capture  Cross  Growth  method 
Section  (cur) 

Ref. 

ET1 

E--0.85 

6.5x10“" 

Bulk 

23 

ET2 

0.30 

2.5xl0“15 

fl 

It 

ESI 

0.83 

1.0X10"13 

Bulk 

24 

Ell 

0.43 

7.3x10“*® 

VPE 

25 

El  2 

0.19 

i.ixur« 

VPE 

It 

El  3 

0.18 

2.2xl0"14 

VPE 

It 

EB1 

0.86 

3.5xl0“*4 

Cr-doped  LPE 

26 

EB2 

0.83 

2.2xl0“13 

as-grown  VPE 

It 

EB3 

0.90 

3.0x10“** 

Electr.irradi . 

27 

EB4 

0.71 

8.3xl0“13 

It  It 

It 

EB5 

0.48 

2.6xl0“*3 

MBE 

28 

EB6 

0.41 

2.6xl0“13 

Electr .  i rradi . 

27 

EB7 

0.30 

1.7xl0“*4 

MBE 

28 

EB8 

0.19 

1.5x10“ 14 

MBE 

28 

EB9 

0.18 

- 

Electr.  i rradi. 

27 

EB10 

0.12 

- 

It  It 

It 

ELI 

0.78 

1.0Xl0“14 

Cr-doped  bulk 

29 

EL2  (A) 

0.83 

1.2xl0“*3 

VPE 

30 

EL3  (B) 

0.58 

II 

II 

II 

EL4 

0.51 

1.0xl0“|2 

MBE 

EL5  (C) 

0.42 

1.2xl0“*3 

VPE 

30 

EL6 

0.35 

1.5xl0“*3 

Bulk 

29 

EL7 

0.30 

7.2x10“*® 

MBE 

EL8  (D) 

0.28 

7.7x10“*® 

VPE 

30 

EL9 

0.23 

6.8x10“*® 

It 

It 

EL10 

0.17 

1.8x10“" 

MBE 

EL11  (F) 

0.17 

3.0x10“*-® 

VPE 

30 

EL12  (A) 

0.78 

4.9x1 0“J‘ 

•I 

II 

EL14 

0.22 

5. 2x10“ *® 

Bulk 

29 

EL15 

0.15 

5. 7x10“ j3 

Electr.  i rradi. 

EL16 

0.37 

4.0xl0“*8 

VPE 

30 

EFl 

0.11 

Proton  i rradi . 

18 

EF2 

0.14 

It  It 

II 

EF3 

0.20 

It  It 

II 

EF4 

0.31 

Electr.irradi . 

21 

EF5 

0.35 

n-LEC 

EF6 

0.52 

Proton  i rradi. 

18 

EF7 

0.60 

n-LEC 

EF8 

0.71 

Electr.  i rradi . 

21 

EF9 

0.76 

LEC 

EF10 

0.83 

LEG,  VPE,  MOCVD 

31 

EF11 

0.90 

Electr .  i rradi . 

21 

» 
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Table  2.3  Hole  traps  in  p-GaAs. 


Level 

Activation 

Capture  Cross 

Growth  method 

Ref. 

energy  (eV) 

Section  (an2) 

HT1 

Ey+0.44 

1.2xl0“14 

VPE 

23 

HS1 

0.58 

2.0x10“^ 

LPE 

24 

HS2 

0.64 

4.1X10T}-® 

II 

It 

HS3 

0.44 

4.8xl0“18 

II 

II 

HBl 

0.78 

5.2xl0"l® 

Cr-doped  LPE 

26 

HB2  (B) 

0.71 

1.2x10“" 

as-grown  LPE 

II 

HB3 

0.52 

3.4x10-}® 

Fe-doped  LPE 

II 

HB4 

0.44 

3.4xl0-~4 

Cu-doped  LPE 

•I 

HB5  (A) 

0.40 

2.2x10-1® 

as-grown  LPE 

II 

HB6 

0.29 

2.0X10-14 

Electr .  j  rradi . 

27 

HL1 

HL2 

0.94 

0.73 

3.7xl0“l4 

1.9xl0“r4 

3.0x10-1® 

3.0xl0“l® 

9. 0x10“  J-4 
5.6xl0’r4 

Cr-doped  VPE 
as-grown  LPE 

32 

HL3 

0.59 

Fe-doped  VPE 

33 

HL4 

0.42 

Cu-doped  VPE 

II 

HL5 

0.41 

as-grown  LPE 

32 

HL6 

0.32 

VPE 

33 

HL7 

0.35 

6.4x10-1® 

MBE 

It 

HL8 

0.52 

3.5x10-1® 

MBE 

II 

HL9 

0.69 

1.1x10-" 

VPE 

II 

HL10 

0.83 

1.7x10-1® 

VPE 

II 

HF0 

0.082 

HF1 

0.13 

HF2 

0.17 

HF3 

0.20 

HF4 

0.29 

HF5 

0.35 

HF6 

0.40 

HF7 

0.44 

HF8 

0.52 

HF9 

0.57 

LEC  (p)  17 

Electr .  i rrad i .  21 

Proton  irradi.  18 

••  II  II 

Electr.  irradi.  2 

II  It  II 

II  II  II 

Proton  irradi.  18 

II 

II 


III.  EXPHUmiEAL  DETAILS 


Experimental  tools  employed  in  this  study  include  cur  rent- voltage  (I- 
V)  measurement,  a.c.admi  ttance  measurement,  capaci  tance- voltage  (C-V) 
measurement,  thermally  stimulated  capacitance  (TSCAP)  measurement,  and  deep 
level  transient  spectroscopy  (DLTS)  measurement.  From  these  measurements, 
one  can  determine  the  diode  characteristics  and  the  defect  parameters  such 
as  defect  energy  levels,  defect  densities,  and  capture  cross  sections. 
Each  of  these  measurement  technique  is  described  as  follows: 

3.1.  Current-Voltage  (I-V)  Measurement: 

Measurement  of  the  current- voltage  (I-V)  characteristics  under  forward 
bias  condition  can  yield  useful  information  concerning  the  conduction 
mechanisms,  recombination  processes  in  the  space  charge  region  of  a  p-n 
junction  diode  or  a  Schottky  barrier  diode.  For  a  good  p-n  diode  with  no 
surface  leakage,  the  total  current  is  composed  of  the  diffusion  current  in 
the  quasi -neutral  region  (QNR)  and  the  recombination  current  in  the  space 
charge  region  (SCR).  When  bulk  diffusion  current  component  dominates,  the 
current  expression  is  given  by:^34^ 

If  -  Id  (exp(Va/VT)  -  1]  (3.1) 

where  Vg  is  the  applied  voltage;  VT  =  kT/q,  and  ld  is  the  magnitude  of 
the  saturation  diffusion  current. 

Id  =  qnj2  A  ((VL^)  +  (Dp/L^)  ]  (3.2) 

where  n^  is  the  intrinsic  carrier  density;  A  is  the  diode  area;  Dn  (Dp)  is 
the  diffusion  constant  for  electrons  (holes);  Ln  (Lp)  is  the  diffusion 
length  for  electrons  (holes),  and  Na  (Nd)  is  the  acceptor  (donor)  density. 
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If  bulk  generation-recombination  current  dominates,  then  the  current 
expression  is  given  by:^4l 


Ir  =  Irg  exp(va/2VT) 


(3.3) 


where  lrg  is  the  magnitude  of  the  generation-recombination  current. 


Irg  *  <r>iW2\ 


(3.4) 


In  Eq.(3.4),  W  denotes  the  depletion  layer  width;  7^  =  (T^Tp)1^  is 
the  effective  carrier  lifetime  in  the  space  charge  region;  ^51  Tn  (7£>) 
is  the  lifetime  of  electrons  (holes),  defined  by: 


%  *  1  /  <Nt  «*n  <vth» 


(3.5) 


where  Nfc  is  the  trap  density;  (  S'p)  is  the  capture  cross  section  for 
electrons  (holes);  <vfch>  is  the  average  thermal  velocity. 

The  total  current  can  be  expressed  by: 


It  -  If  +  Ir  -  IQ  exp(Va/nVT) 


(3.6) 


where  IQ  is  the  saturation  current;  n  is  the  ideality  factor  of  a  p-n  diode 
which  is  usually  used  to  identify  the  dominant  current  component  in  a  p-n 
diode.  Inspection  of  Eqs.(3.1)  and  (3.3)  shows  that  the  bulk  diffusion 
current  depends  more  strongly  on  temperature  than  the  recombination  current 
in  the  junction  space  charge  region  (SCR).  Since  the  recombination  current 
in  the  junction  SCR  is  inversely  proportional  to  the  effective  carrier 
lifetimes  (and  hence  directly  related  to  the  defect  density  in  the 
transition  region),  measurements  of  I-V  characteristics  under  forward  bias 


conditions  would  allow  us  to  determine  the  effective  lifetimes  in  GaAs. 


3.2.  A.  C.  Admittance  Measument. 

A  useful  experimental  tool  for  evaluating  shunt  (Rp)  and  series 
resistance  (Rg)  across  a  p-n  junction  diode  is  by  using  a.c.admi ttance 
measurement  technique. ^*>1  a  p-n  junction  diode  can  be  represented  by  a 
three  element  device  with  shunt  resistance,  junction  capacitance,  and 
series  resistance.  Measurements  of  the  admittance  as  a  function  of 
frequencies  (e.g.  from  110  kHZ  to  700  MHZ)  will  enable  us  to  determine 
values  of  R_,  R_,  and  C  in  a  p-n  diode.  The  impedence  of  the  p-n  junction 
diode  is  given  by: 

1 

Z(w)  =  Rs  +  {Rp(l/jwC)/[Rp  +  (i/jwC)]} 

-  Rjjtw)  +  jXD(w)  (3.7)  _  , 

where 

R^w)  =  Rs  +  Rp/  ( l+w2C2Rp2 )  (3.8) 

XD(w)  =  RpwC/(l4w2C2Rp2)  (3.9) 

The  admittance  is  the  inverse  of  the  impedence  in  a  p-n  junction  diode, 
which  is  given  by: 

Y  (w)  =  1/Z  (w)  =  (Rd(w)/(Rd2(w)  +  Xd2(w)]}  ‘  1 

'  -1 

-  j(XD(w)/[R02(w)  +  Xd2(w)]}  (3.10)  ; 

In  the  low  frequency  and  high  frequency  limits: 

w  ->  0,  Re  Y (w)  =  1  /  (Rp  +  Rs]  (3.11) 

w  ->  oo  ,  Re  Y (w)  =  1  /  Rs  (3.12) 

Thus,  from  a  plot  of  the  Im  Y(w)  vs.  Re  Y(w),  the  series  resistance  and 


shunt  resistance  of  the  diode  can  be  determined.  The  a.c.admi  ttance 


technique  can  he  very  useful  in  evaluating  the  shunt  leakage  problems  in  a 

p-n  diode  and  allows  accurate  determination  of  the  Re,  R_,  and  C  over  a 

s  p 

wide  range  of  frequencies. 

3.3.  Capacitance-VOltage  (C-V)  Measuraent: 

The  capaci tance- voltage  (C-V)  measurement  can  be  used  to  determine  the 
background  doping  concentration  in  the  n-  or  p-GaAs  epitaxial  layers  using 
a  Schottky  barrier  structure  or  a  one-sided  abrupt  p+-n  (or  n+-p)  junction. 
The  depletion  capacitance  across  the  Schottky  barrier  diode  is  given  by: 

C(Vr)  *  6gA/W  *  A  {qtfsNd/[2(4j  +  Vr  -  kT/q)] }1/2  (3.13) 

where  fig  is  the  dielectric  constant  of  GaAs;  is  the  built-in  potential; 
and  Vr  is  the  applied  reverse  voltage.  Eq.(3.13)  shows  that  the  depletion 
capacitance  of  a  Schottky  diode  is  proportional  to  the  square  root  of 
dopant  concentration  and  inversely  proportional  to  the  square  root  of  the 
applied  voltage.  If  the  inverse  of  the  capacitance  square  (c-2)  is  plotted 
as  a  function  of  the  reverse  bias  voltage  Vr,  then,  the  background 
concentration,  can  be  calculated  from  the  slope  of  C-2  vs  Vr  using  the 
following  expression: 

C"2(Vr)  =  (2  /  (q  6gA2Np)  J  (^  +  Vr)  (3.14) 

The  intercept  of  C“2  vs.  Vr  plot  in  the  voltage  axis  yields  values  of  Oj 
which  is  related  to  the  barrier  height  of  a  Schottky  diode  by: 

4sn  =  $i  +  vn  +  kT/<J  -  (3.15) 

where 

vn  =  Ej,  -  (kT/q)  ln(Nd/Nc)  (3.16) 
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and  is  the  image  lowering  potential  of  a  Schottky  diode;  Nc  is  the 
effective  density  of  states  in  the  conduction  hand. 


3.4.  Ttaemally  Stimulated  Capacitance  (T9CAP)  Measurement: 

Another  interesting  experiment,  which  is  known  as  the  thermally 
stimulated  capaci tance  (TSCAP)  measurement  technique  137,38]  wjii  be 
described  in  this  section.  The  TSCAP  experiment  is  carried  out  by  first 
reverse  biasing  a  p-n  diode  or  a  Schottky  diode,  and  then  the  diode  is 
cooled  down  to  liquid  nitrogen  temperature  (77°k).  After  temperature 
reaches  77°K,  the  diode  is  momentarily  zero  biased  to  fill  the  majority 
carrier  traps  and  returned  to  reverse  bias  condition,  and  the  temperature 
is  then  raised  from  77  to  400  °K.  The  thermal  scan  of  capacitance  vs 
temperature  plot  is  then  taken  by  using  an  X-Y  recorder.  A  capacitance  step 
is  observed  from  the  C  vs  T  plot  if  majority  or  minority  carrier  emission 
is  taking  place  in  a  trap  level  within  a  certain  temperature  range.  The 
amplitude  of  this  capacitance  step  is  directly  proportional  to  the  trap 
density.  The  trap  density  for  an  n-GaAs  can  be  calculated  from  the 
following  expression: 

Nt  =  Nd(2AC  /  CQ)  (3.17) 

where  CQ  is  the  depletion  layer  capacitance  and  AC  is  the  capacitance 
change  due  to  the  majority  or  minority  carrier  emission.  Thus,  knowing  Nd 
(or  Ng)  and  CQ  at  the  temperature  where  the  capacitance  step  was  observed, 
the  trap  density  can  be  calculated  from  Eq.(3.17).  Note  that  Eq.(3.17)  is 
valid  only  for  the  case  when  Nfc  is  less  than  0.1  Nd.  For  the  case  of  large 
trap  density  with  Nfc  >  0.1  Nd,  a  more  exact  expression  should  be  used 


instead 


3.5.  DLTS  Measurement: 


The  Deep-Level-Transient-Spectroscopy  (DLTS)  experiment  is  a  high 
frequency  (20  MHZ)  transient  capacitance  technique,  which  was  introduced 
first  by  Lang  in  1974  [7,14].  T^e  DLTg  scan  displays  the  spectrum  of  deep 
level  traps  in  the  forbidden  gap  of  a  semiconductor  as  positive  or  negative 
peaks  on  a  flat  baseline  as  a  function  of  temperature.  Although  this  kind 
of  measurement  is  time  consuming,  it  offers  several  advantages  such  as 
sensitive,  easy  to  analyze  and  capable  of  measuring  the  traps  over  a  wide 
range  of  depth  in  the  forbidden  gap.  By  properly  changing  the  experimental 
conditions,  one  can  measure  the  following  parameters: 

-  Minority  and  majority  carrier  traps. 

-  Activation  energy  of  each  defect  level. 

-  Defect  concentration  which  is  directly  proportional  to 
the  peak  height. 

-  Defect  concentration  profile. 

-  Electron  and  hole  capture  cross  sections. 

3.5.1.  Principles  of  title  DLTS  Technique. 

The  capacitance  transient  is  associated  with  the  return  to  thermal 
equilibrium  of  the  carrier  occupancy  in  a  trap  level  following  an  initial 
nonequilibrium  condition.  The  polarity  of  the  DLTS  peak  depends  on  the 
capacitance  change  after  trapping  the  minority  or  majority  carriers. 
Because  an  increase  in  trapped  minority  carriers  in  the  SCR  would  result  in 
an  increase  in  the  junction  capacitance,  the  minority  carriers  trapping 
will  produce  a  positive  polarity  peak,  and  vice  verse.  For  example,  in  a 
p+n  junction  diode,  the  SCR  extends  mainly  into  the  n-region,  and  the  local 
charges  are  due  to  positively  charged  ionized  donors.  If  a  forward  bias  is 


allied,  the  minority  carriers  (holes)  will  be  injected  into  this  region. 
Once  the  holes  are  trapped  in  a  defect  level  the  net  positive  charges  in 
such  region  will  increase.  This  results  in  a  narrower  SCR  width  which 
implies  a  positive  capacitance  change.  Thus,  the  DLTS  signal  will  have  a 
positive  peak.  Similarly,  if  the  majority  carriers  are  injected  into  this 
region  and  captured  by  the  majority  carrier  traps,  which  reduces  the  local 
charges,  the  SCR  width  will  be  wider,  implying  a  decrease  of  the  junction 
capacitance.  Therefore,  the  majority  carrier  trapping  will  result  in  a 
negative  DLTS  peak.  The  same  argument  can  be  applied  to  the  n+p  junction 
diodes.  All  of  the  samples  used  in  the  DLTS  measurements  are  p+n  diodes 
so  that  a  positive  peak  represents  the  hole  trap  and  a  negative  peak 
represents  the  electron  trap. 


3.5.2.  Defect  Concentration 

The  defect  concentration  is  directly  proportional  to  the  peak  height 
as  described  before,  and  the  peak  height  is  proportional  to  the  capacitance 
change  C(0).  Therefore,  the  defect  concentration  Nt  is  proportional  to 
C(0)  which  can  be  derived  as  follows:  Let  C(t)  be  the  capacitance 
transient t39'40!  which  is  proportional  to  electrons  (hole)  emitted  to  the 
conduction  (valence)  band,  then, 

C(t)  =  A{q  fis  [Nd  -  Ntexp(-t/r)l  /  [2($j  +  Vr  +  kT/q)]}1/2 
=  CQ  {1  -  [Ntexp(-t/r)  /  Nd]}1/2  (3.18) 

where  t  is  time;  is  the  carrier  emission  time  constant;  CQ  =  C(Vr),  as 
shown  in  Eq.(3.13),  is  the  junction  capacitance  at  the  quiescent  reverse 
bias  condition.  Using  binomial  expansion  and  the  condition  that  Nt/tod  « 
1,  Eq.  (3.18)  reduces  to  a  simple  form  as: 

C(t)  =  CQ  (1  -  Ntexp(-t/T)/2Nd]  (3.19) 
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Eq. (3.19)  can  be  rewritten  as: 


Ntexp(-t/r)  =  (2iC(t)  /  CD)  Nd  (3.20) 

where4C(t)=CQ-C(t).  From  DLTS  measurement,  C(0)  can  be  determined.  The 
junction  capacitance  CQ  and  the  background  concentration  can  be  obtained 
from  C-V  measurements.  Thus,  the  defect  concentration  Nfc  can  be  calculated 
easily  by  using  Eq.(3.20)  at  t=0. 

3.5.3.  Activation  Energy  of  the  Defect  Level 

The  decay  time  constant  in  the  capacitance  transient  during  the  DLTS 
experiment  is  associated  with  a  specific  time  constant  which  is  equal  to 
the  reciprocal  of  the  emission  rate.  For  an  electron  trap,  the  emission 
rate  "en"  is  functions  of  temperature,  capture  coefficient  and  activation 
energy,  and  can  be  expressed  by^'^1 

en  *  (°”n<vth>Nc  /  exp((Ec  -  Et)/kT]  (3.21) 

where  Et  is  the  activation  energy  of  the  trap;  g  is  the  degeneracy 
factor;  <Sn  is  the  electron  capture  cross  section  which  is  temperature 
dependent t 43 1 ,  and  is  given  by 

<Tn  =*<3^  expf-AEt/kT)  (3.22) 

where <$*00  is  the  capture  cross  section  at  very  high  temperature;  AE^  is 
the  activation  energy  of  the  capture  cross  section  for  the  trap.  en 
can  be  written  as: 

en  =  BT2  exp{[Ec  -  (Efc  +  4Eb)]  /  kT} 

=  BT2  exp  [(Ec  -  V  /  kTl  <3-23) 
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where  B  is  the  proportionality  constant,  and  is  independent  of  temperature. 
From  this  relation,  en  increases  with  increasing  temperature.  The 
capacitance  transient  is  rearranged  from  Eq.(3.19)  as: 

C(t)  =  C0(Nt/2Nd)  exEK-t/r) 

=  4C(0)  exp(-t/r)  (3.24) 

where  ®n_1/  is  the  reciprocal  emission  time  constant. 

The  procedures  for  determining  the  activation  energy  of  a  defect 
level  in  a  semiconductor  is  described  as  follows.  We  first  set  t^  and  t2 
in  a  dual-gated  integrator  boxcar,  then  we  can  write: 

C(tL)  =  C(0)  expl-tj/f)  (3.25) 

C(t2)  =  C(0)  exp(-t2/r)  (3.26) 

The  DLTS  scan  along  the  temperature  axis  is  obtained  by  taking  the 
difference  of  Eq.(3.25)  and  (3.26),  which  yields: 

S(T)  =  C(0)  [expt-tj/r)  -exp(-t2/T)]  (3.27) 

The  maximum  emission  rate,  Tmax*1,  is  obtained  by  di  fferentiating  SfZf 
with  respect  to '2',  and  sets  dS(T)/d  T  »  0,  which  yields: 

Tmax  =  <fcl  -  0.28) 

Under  this  condition,  S(tf  reaches  its  maximum  value  at  a  specific 
temperature.  The  emission  rate  is  given  by  en  =  V"2Tmax  for  each  t^  and 
t2  setting.  By  changing  t^  and  t2  several  times,  a  set  of  temperatures 
correspond  to  this  set  of  mgx  (or  emission  rate  en)  can  be  obtained  as 
shown  in  Fig.3.1.  The  activation  energy  of  the  trap  can  be  calculated  from 
the  slope  of  the  Arrhenius  plot. 
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IV.  DBTBMINAXIGN  OF  POTBfTIAL  WHi  OF  DfflP-LEVH.  TRAPS  IN  GAAS 

Deep-Level  Transient  Spectroscopy  (DLTS)  signals  are  often  analyzed 
assuming  that  the  capacitance  transient  is  exponential.  However,  most  of 
the  DLTS  signals  are  nonexponential  transient^^'^ .  Nonexponential 
capacitance  transients  may  occur  from:  (1)  electric  field  dependent 
emission  rate  of  the  trapped  charge  ^  44,4S^ ,  (2)  multi  exponential  decay  due 
to  several  trap  levels  with  similar  emission  rates,  and  (3)  trap  density 
of  the  same  magnitude  as  that  of  the  shallow  dopant  density^48'4^.  This 
chapter  will  deal  with  cases  (1)  and  (2). 

To  determine  the  nonexponential  capacitance  transient  due  to  electric 
field  dependent  emission  rate  of  the  trapped  charge,  Makran-Efrid ^ 48,49 1 
and  Wang  et  ai [50,51]  measured  the  field  enhancement  emission  rate. 
In  this  chapter,  we  present  the  theoretical  analysis  of  nonexponential 
capacitance  transients  due  to  the  electric  field  dependence  emission  rate 
of  trapped  charge  and  multi  exponential  decay.  Emission  rate  of  trapped 
charge  carriers  is  enhanced  by  the  Poole-Frenkel  effect  and  phonon-assisted 
tunneling  effect  in  the  presence  of  an  applied  electric  field.  Since  the 
electric  field  varies  with  position  within  the  depletion  region  of  a 
reverse  biased  p-n  junction,  the  emission  rate  is  not  constant  within  the 
same  region.  The  DLTS  response  due  to  the  nonuni  formed  emission  rates  can 
in  general  be  expressed  by  S ('Q  =  2  exp(-enjt^)  -  £  exp(-enjt2).  Based 
on  the  nonexponential  capacitance  transient,  theoretical  calculations  of 
DLTS  response  for  the  deep- level  traps  in  GaAs  are  made  using  five 
different  potential  wells;  namely,  the  Coulombic  potential  well  which  has  a 
positive  charge  state  for  the  empty  state,  the  Dirac  well,  square  well, 
polarization  well  as  well  as  dipole  well  which  all  have  a  neutral  charge 
state  from  different  physical  origins.  A  comparison  of  the  theoretical 
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calculations  of  nonexponential  DLTS  response  with  the  DLTS  data  for  each 
trap  level  would  allow  us  to  determine  the  potential  well  for  the  trap 
involved.  This  model  was  applied  to  the  calculations  of  potential  well  for 
the  EL2  level  and  other  electron  traps  in  GaAs.  We  found  that  the  EL2 
electron  trap  was  due  to  Coulombic  potential  well  with  a  double  charge 
state  (i.e.,  AsGa++).  Details  of  the  results  will  be  discussed  in  this 
chapter  6.. 

4.1  Capture  and  mission  process  at  a  Deep-Level  Trap. 

in  the  study  of  electric  field  dependence  of  emission  rate,  it  is 
important  to  know  the  charge  state  of  a  trap  so  that  the  type  of  potential 
well  for  such  a  trap  can  be  determined.  The  charge  state  of  a  deep- level 
trap  may  be  either  positively  charged,  neutral,  or  negatively  charged ^ 10 J 
(i.e.,  Nt+,  Nt°  or  Nt").  For  example,  if  an  electron  is  captured  by  a 
positive  donor  trap,  then  the  kinetic  equation  for  the  capture  process 
requires  that: 

Nt+  +  e  =  Nt°.  (4.1) 

In  this  case,  the  capture  process  is  Coulombic  attractive.  Similarly,  the 
emission  process  for  the  same  trap  can  be  written  as 

Nt°  *  Nt+  e,  (4.2) 

which  again  is  an  attractive  process.  This  attractive  capture- 
emission  process  is  the  feature  of  a  Coulombic  potential  well.  In 
GaAs,  for  example,  the  deep-level  defects  such  as  arsenic  vacancy 
1*1(10  and  arsenic  antisite^l  (AsGa‘H‘)  defects  have  a  Coulombic 
potential  well. 
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For  an  empty  neutral  trap  with  Nt°,  the  capture  process  is  given  by 

Nt°  +  e  =  Nt“f  (4.3) 

,  and  the  emission  process  for  the  same  trap  is  given  by 

Nt“  =  Nt°  +  e.  (4.4) 

The  potential  wells  for  a  deep-level  neutral  trap  in  a  semiconductor  may 

include  Dirac  delta  potential  well^53^,  square  potential  well,  polarization 

potential  well  and  dipole  potential  well[54].  The  Dirac  delta  potential 

well  does  not  show  Poole-Frenkel  effect.  Neutral  arsenic  vacancy  (VAs°) 

defect  in  a  GaAs  crystal  behaves  as  a  Dirac  delta  potential  well  in  the 

capture-emission  process.  For  the  square  potential  well,  its  potential 

well  has  a  depth  of  VQ  for  r  <  rQ  and  zero  for  r  >  rQ.  Deep- level  defect 

such  as  neutral  arsenic  interstitial  (ASj°)  has  this  characteristics.  The 

4 

potential  of  a  polarization  potential  well  has  the  form  V(r)  =  -  A-^  /  r  , 
where  Aj^  =  e2*  j  /8^>€s2»  is  the  polarizability  of  a  neutral  impurity 
atom,  and  6S  is  the  relative  dielectric  constant  of  the  host  crystal. 

t 55]  yged  a  polarization  potential  well  to  model  capture  of  electrons 
and  holes  by  neutral  impurities.  Tasch  and  Sah^44^  used  this  potential 
well  to  model  the  observed  field  dependence  of  Au  in  silicon.  Neutral 
impurity  trap  such  as  Al^0  has  the  properties  of  a  polarization  potential 
well.  A  complex  of  two  oppositely  charged  centers,  each  with  charge 
magnitude  of  Zq,  will  form  a  dipole  potential  well.  The  potential  is 
represented  by  a  dipole  potential,  V(r)  *  qr • P/4fiQfirs2,  with  a  dipole 
moment  given  by  P  =  Z  q  RQ  ,  where  RQ  is  the  distance  separating  the  two 
charge  centers.  The  antisite  pair  defect  (AsGa++GaAs  )  behaves  as  a 
dipole  potential  well  in  the  capture-emission  process. 
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For  the  process  Nt“  +  e  =  Nt  ,  the  trapping  process  is  repulsive  and 
the  trap  in  Nt“  state  has  a  Coulombic  energy  harrier  that  the  electron  must 
surmount  thermally  or  tunnel  through.  Thus  one  may  expect  i ts  capture 
cross  section  to  be  very  small  (~10  cm*  or  less)  and  with  a  strong 
temperature  dependent  property. 

In  DLTS  technique ^»14]  f  electrons  captured  by  a  trap  level  are 
occured  in  the  quasi -neutral  region  (QNR)  of  a  p-n  diode  which  has  zero 
electric  field.  On  the  other  hand,  electrons  emission  from  a  trap  level 
usually  take  place  in  the  space-charge  region  (SCR)  of  a  p-n  diode  in  which 
high  electric  field  prevails  in  this  region.  Fig. 4.1  shows  that  the 
capture  process  in  QNR  and  emission  process  in  SCR.  The  effect  of  electric 
field  on  the  emission  process  is  to  enhance  the  emission  rate  by  either 
Poole-Frenkel  effect  or  phonon-assisted  tunneling  effect  to  be  discussed  in 
next  section.  The  emission  rate  enhanced  by  the  electric  field  will  result 
in  a  nonexponential  DLTS  response. 

4.2  Overview  of  the  Theory. 

The  emission  rate  is  based  on  the  detailed  balence  expressi on ^ 41,42 ^ , 
and  under  zero  electric  field  condition,  Eq.(3.23)  is  rearranged  as: 

eno  *  %  <»exP(-  EtjAT)  (4.5) 

where  Eti  is  the  ionization  energy  of  the  trap  level,  and  en  contains  the 
matrix  elements  of  the  transition.  Fig.4.2  shows  a  trap  level  under  high 
electric  field.  There  are  three  basic  mechani sms ^4^  which  affect  the 
emission  rate  enhancement  under  high  electric  field  conditions. 

(1)  Poole-Frenkel  effect,  in  which  electrons  climb  over  a 

barrier  lowered  by  the  presence  of  applied  electric  field. 

(2)  Phonon-assisted  tunneling  effect,  in  which  electrons  absorb 


24 


25 


thermal  energy  from  the  lattice,  and  then  tunnel  through  the 
harrier  at  a  higher  energy. 

(3)  Pure  quantum  mechanical  tunneling  effect.  Since  pure  tunneling 
becomes  important  only  at  very  high  electric  fields ^*>1  (j.e.  p 

>  107  V/cm) ,  we  will  not  consider  this  effect  in  our  DLTS 
analysis. 


4.2.1  Poole-Frenkel  effect. 

The  Poole-Frenkel  effect  is  associated  with  the  lowering  of  potential 
barrier  of  a  deep-level  trap  in  a  bulk  semiconductor.  Both  the  donor  and 
acceptor  traps  wi  11  show  the  Poole-Frenkel  effect.  The  enhancement  of 
emission  rate  from  a  Coulombic  potential  well  due  to  Poole-Frenkel  effect 
was  first  done  in  a  one-dimensional  model  by  Frenkel ^57J,  and  later 
extended  to  three-dimensional  case  independently  by  Hartke^58^  and 
Jonscher^59^.  A  similar  calculation  for  the  polarization  well  and  dipole 
well  was  done  by  Martin  et  al^54^. 

In  this  section  we  will  consider  the  general  case  for  calculating  the 
emission  rate  enhancement  due  to  Poole-Frenkel  effect.  The  potential  for  a 
deep-level  trap  in  an  electric  field,  F,  can  be  expressed  by: 

VT(r)  *  V(r)  -  q  F  r  cos(0)  (4.6) 

where  V(r)  is  the  potential  of  a  deep-level  trap.  For  0  <  &  <  k/2,  the 
minimum  potential  was  found  by  setting  dVT(r)/dr  =  0  at  r  =  rmax.  The 
lowering  of  potential  barrier  due  to  the  presence  of  an  applied  electric 
field  is  found  by  evaluating  VT(r)  at  r  =  rnax.  Thus, 

*EtiW)  -  VT(rmax)  (4.7) 
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Fig. 4. 2  Enhancement  of  electron  emission  rate  by:  (l)  Poole-Frenkel 

effect,  (2)  phonor-assisted  tunneling  effect,  (3)  pure 
quantum  mechanical  tunneling  effect. 
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One  dimensional  result  is  given  by  setting  9  =  0.  The  Poole-Frenkel 
effect  leads  to  a  decrease  in  ionization  energy,  and  the  one  dimensional 
thermal  emission  rate,  enp  in  the  presence  of  an  electric  field  can  be 
written  as 


®ni  =  exp[-(Etj  -  AEfci)  /  kTj 

■  enQexp(  EtiAT)  (4.8) 

where  eno  is  the  thermal  emission  rate  at  zero  electric  field  given  by 
Eq.(4.5).  The  three-dimensional  calculation  requires  an  integration  over 
0  due  to  the  spatial  variation  of  Efcj  (  ).  For  0  <  6  <  7f/2,  the  emission 
rate  is  proportional  to  0  ,  and  for  7T/2  <0<7t,itis  assumed  that  the 
emission  rate  does  not  change  with  electric  field.  Thus, 

en3/eno*(l/4fl)[ /^d$J^sin(e)exp(4Ej  AT)d®  +/^d$jQsin(®)  de]  (4.9) 

where  is  the  three  dimensional  thermal  emission  rate. 

4.2.2  Ptionan-^assi  stood  tunneling  effect. 

There  are  several  different  types  of  potential  wells  which  may  exist 
in  a  deep-level  trap,  depending  on  the  charge  state  of  a  particular  trap. 
The  potential  well  may  be  affected  by  an  external  electric  field.  Fiq. 4.3 
shows  a  potential  well  for  a  trap  level  with  energy  of  Etj  as  a  function  of 
the  radial  coordinate.  A  trapped  electron  can  absorb  a  phonon  and  tunnel 
through  the  barrier.  The  tunneling  process  can  be  treated  by  using  the 
(Wentzel-Kramers-Bri  llouin)  approximation  for  the  potential 
barrier.  The  probability  ratio  (or  the  barrier  penetration  factor,  P) 

for  an  electron  impinging  on  the  harrier  is  given  by 

x2 

P  *  exp[-2  f  k(x)  dx] 

X1 
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(4.10) 


=  exp{-[  (8m*)  ^2/)(]  J  2[V(x)-E]1//2dx} 

X1 

where  k(x)  =  ( (2m*)  1//2/>i]  x  [V(x)-E]^2  is  positive  when  V(x)  >  E.  V(x)  is 
the  potential  well  for  the  trap  level.  E  is  the  electron  energy,  x-^  and 
x2  are  the  two  turning  points  for  which  electron  energy  is  equal  to  the 
barrier  energy  of  the  trap  potential  well. 

A  trapped  electron  can  absorb  a  phonon  with  energy  of  exp(-EfchAT)  and 
tunnel  through  the  barrier  at  a  higher  energy.  The  probability  of  this 
composite  event  is 

Pc  =  exp(-Eth/kT)  P  (4.11) 

where  E^  is  the  thermal  excitation  energy  or  phonon  energy;  and  P  is  given 
by  Eq.(4.10).  The  emission  rate  due  to  phonon-assisted  tunneling  effect  is 
given  by  integrating  Eq.(4.11)  over  the  phonon  energy,  Eth. 

Ej  AT 

eht  *  f  pc  dEthAT  (4.12) 

' o 

where  1/kT  js  a  normalization  factor.  The  normalized  emission  rate 
enhancement  is  obtained  by  dividing  Eq.(4.12)  by  the  zero  electric  field 
emission  rate  e^,  which  yields 

EjAT 

®nt/eno  *  exP(EjAT)  (  Pc  dEtl/kT 

o 

E  •  /kT 

*  exp[(Eti-Eth)AT]  (  3  exp{  [—(8m*)  1//2/h] 

Jq 

x  [V(x)-E]dx  dEthAT}  (4.13) 

Eq.(4.13)  is  a  general  expression  for  the  emission  rate  enhancement  due  to 
phonon-assisted  tunneling  effect. 
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The  total  emission  rate  enhancement  is  obtained  from  the  sum  of  Poole- 


Frenkel  effect  and  phonon-assisted  tunneling  effect.  From  Eq.(4.9)  and 
(4.13),  one  obtains: 

enH  =  ^0*^3  +  ®ht>  <4*14) 

Vincent  et  al^3]  and  Martin  et  al^4^  showed  that  for  a  Coulombic  well  and 
other  types  of  potential  wells,  both  the  Poole-Frenkel  effect  and  phonon- 
assisted  tunneling  effect  are  important  over  the  field  range  of  interest 
(106  -  108  V/m  at  300  K).  Our  results  are  in  good  agreement  with  this 
observation,  as  will  be  shown  next. 

4.3  Theoretical  calculations  of  the  emission  rate  fnhanoMcnt 
for  different  potential  wells. 

4.3.1  Couloebic  potential  well. 

A  donor  type  electron  trap,  which  has  a  Coulombic  potential  well,  is 
positively  charged  (e.g„  Nfc+)  when  it  is  empty.  The  capture  and  emission 
process  for  a  such  trap  is  given  by  Eq.(4.1)  and  (4.2),  respectively. 
Assuming  an  electron  is  at  a  distance  ,  r  ,  from  the  trap,  the  attractive 
force  between  a  positive  trap  and  the  electron  is  given  by 

F(r)  =  -q2  /  4*^s2  (4.15) 

The  potential  for  the  trap  is  obtained  by  integrating  Eq.(4.15)  from 
infinity  to  r.  Thus, 

V(r)  «/  F(r)  dr  =  -qV4  *60egr  (4.16) 

When  an  external  field,  F,  is  applied,  the  total  potential  energy  as  a 
function  of  distance  is  given  by 

Vm(r)  *  -q2/4 n  60<rr  -  qFrcos(e) 


(4.17) 


)  defects 


Electron  traps  such  as  As  vacancy  (Vfts+)  and  As  anti  site  (ASQa++ 
have  coulombic  potential  well  for  electrons.  Eq.  (4.17)  can  also  be 
applied  to  an  acceptor  type  hole  trap,  which  also  has  a  Coulombic  potential 
well.  In  this  case,  the  Coulombic  potential  well  is  negatively  charged  when 
it  is  empty.  The  capture  process  for  such  a  hole  trap  is 

Nt~  +  h+  =  Nt°  (4.18) 

and  the  emission  process  is 

Nfc°  =  Nt“  +  h+  (4.19) 

Hole  traps  such  as  Ga  vacancy  (V^-)  and  Ga  anti  site  (Gafts  )  defects 
have  Coulombic  potential  well  for  holes.  Compare  to  other  types  of 
potential  wells,  Coulombic  potential  well  shows  a  very  large  Poole- 
Frenkel  effect. 

4. 3. 1.1  Poole-Frenkel  effect: 

The  Poole-Frenkel  lowering  potential  Etj  and  the  location  of  lowering 
rmax  are  9*ven  the  condition  dVT(r)/dr  =  0,  or 

rmax  =  {q/[4»60«sF  cos(e)]}1/2 

Etj  =  q  [  q  F  cos(9)/*  £Q€r] (4.20) 

The  Poole-Frenkel  effect  leads  to  a  decrease  of  ionization  energy  Et ^ , 
and  the  one-dimensional  thermal  emission  rate  is  given  by  Eq.  (4.8). 

For  three-dimensional  case,  the  emission  rate  of  a  trapped  electron 
may  be  obtained  by  assuming  that  emission  rate  is  field  dependent  for  0  <9<7T 
/2  ,  and  is  independent  of  electric  field  for  */2  <  &  <7T.  Integrating 
Eq. (4.9)  yields: 
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(4.21) 


en3/eno  =  d/R2)  teR(r-1)  +  1J  +  1/2 

where 

R  =  q  (qF/7Ti0€s)1/2AT  =  ^EtiAT 


(4.22) 


4.3.1.2  Phonon-assisted  txmeling  effect: 

The  phonon-assisted  tunneling  effect  may  be  evaluated  by  substituting 
Eq.(4.17)  to  Eq.(4.13).  Vincent  et  al^l  has  obtained  a  closed  form  for 
the  phonon-assisted  tunneling  effect  for  a  Coulombic  potential  well,  which 
reads: 


(Etj-4Eti)AT 

ent/eno=  j  exp{2-z3//2[4(2m*)1/2(kT)3//2/3qbF] 

o 

x  [l-(aEti/zkT)5/3]}dz  (4.23) 

The  total  emission  rate  consists  of  Poole-Frenkel  and  phonon-assisted 
tunneling  effects,  which  is  expressed  by: 

Sho  =  So  <S3  +  St>  <4-24) 

Fig.  4.4  shows  the  normalized  enhanced  emission  rate  vs  electric  field  for 
the  Coulombic  potential  well  for  the  Ec  -  0.35  eV  electron  trap  oberved  in 
GaAs.  It  is  noted  that  Poole-Frenkel  effect  dominates  when  electric  field 
is  lower  than  2  X  104  V/cm,  while  phonon-assisted  tunneling  effect  becomes 
important  for  electric  field  higher  than  2  x  104  V/cm. 

4.3.2  Dirac  Delta  Potential  well 

The  potential  for  a  Dirac  delta  well  is  given  by: 


i 


Jt. 
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Normalizion  Emission  Rate 


for  r  =  0 


(4.25) 


V(r)  =  -VQ 

=  0  for  |  r  |  >  0 

In  the  presence  of  an  applied  electric  field,  the  potential  is  given  by 
V(r)  =  -VQ  -  qFr  cos(£>  for  r  =  0 

(4.26) 

=  -  qFr  cos(C)  for  |r|  >0 

Neutral  vacancy  trap  such  as  neutral  arsenic  vacancy  (V^0) ,  may  possess  a 
Dirac  potential  well.  For  this  potential  well,  the  Poole-Frenkel  lowering 
potential  is  equal  to  zero,  and  only  the  phonon-assisted  tunneling  effect 
is  considered.  For  the  case  of  a  one-dimensional  calculation,  cos(S)  is 
equal  to  1,  and  the  emission  rate  enhancement  is  obtained  by  substituting 
Eq.  (4.26)  into  Eq.  (4.13),  which  yields: 

£  •  /kT 

ent/eno  a  {exp((Etj-Efch)/kT] }  (  t 1  {expMSm*)1/2/^] 

o 

f /  1  (-qFx  +  Etj  -  Eth)1/2  dx]  dEth/kT}  (4.27) 
x2 

Eq. (4.27)  gives  an  analytical  form  for  the  transparency  of  triangular 
barrier  of  a  Dirac  well.  For  the  triangular  well,  x^  =  0.  and  x2  *  (Eti~ 
Eth)AlF.  Integrating  Eq.(4.27)  over  Efch  and  let  z  =  (Etj  -  Eth)/kT,  we 
obtain  the  enhanced  emission  rate  for  the  phonon-assisted  tunneling  effect 

E  ’/kT 

e^/e^  *Jtl  exp{z  -  z3y/2  [4  (2m*) 1//2  (kT)  3//2/3q)<F] }  dz  (4.28) 
o 

The  total  emission  rate  due  to  Dirac  well  can  be  expressed  by: 

enhr  *  eno<1  +  ent>  <4-29) 


0 
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A33  Square  Potential  Hell 


The  same  treatment  can  he  applied  to  the  square  well.  The  potential 
of  a  square  well  in  the  presence  of  an  applied  electric  field  is 

V(r)  =  V  -  qFr  cos(0)  for  r  <  r_ 

(4.30) 

=  -  qFr  cos($)  for  r  >  rQ 

where  rQ  is  the  radius  of  the  trap  potential  well.  Neutral  interstitial 
trap  such  as  neutral  arsenic  interstitial  (As^°) ,  may  have  such  a  potential 
well.  The  harrier  lowering  due  to  Poole-Frenkel  effect  is  equal  to  qFrQ. 
The  enhanced  emission  rate  for  a  one-dimensional  and  a  three-dimensional 
Poole-Frenkel  effect  is  given  respectively  by: 

’  exp(R’ 

en3/eno  *  d/2  R)  [exp(R)  -  1]  +  1/2 
where  R  *  qFrc/kT. 

For  the  enhanced  emission  rate  due  to  phonon-assisted  tunneling 
effect,  Eq.(4.27)  is  still  valid.  In  a  square  well,  =  rQ  and  x2  =  (Efcj 
-  Eth)/kT.  The  range  of  phonon  energy,  Efch  is  from  0  to  (Et j—  Etj)AT 
Thus,  we  obtain  an  expression  for  the  normalized  enhanced  emission  rate  due 
to  phonon-assisted  tunneling  effect: 

/®ti“  A^ti  )/^T  7/9  *1/9  7/9 

ent/eno  *J  exp{  z-  z3/2[4(2m  )  1/2(kT)  3/2/3qKF] 

yo 

x  [1  -(qFrQ/zkT)3/2]}dz  (4.32) 

The  total  emission  rate  for  the  square  well  is  equal  to  the  sum  of 
Eq.(4.31)  and  Eq.(4.32),  which  is  given  by: 
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enhs  =  eno  ^3  +  ®nt^ 


(4.35) 


4.3.4  Polarisation  potential  well 

An  electron  with  charge,  e,  at  a  distance  r  from  a  trap  center  with 
polarizability, oi  j,  will  induce  a  dipole  moment,  p  =o<je/esr2.  This  dipole 
will  produce  an  attractive  force  on  the  charge  of  pe/2fiQ^r°  *  <XjeV(2 
fcg2r5).  Thus,  the  attactive  potential  is  given  by: 

V(r)  =  -A]/  r4  (4.34) 

where  Aj^  =  fcQ€rs2) 

The  polarizability  of  the  atom  can  be  expressed  by 

<*i/<*H  -  (mo/mJtEj/E^)2  (4.36) 

where  Eti  is  the  ionization  energy  of  the  trap;  EH  *  13.6  eV  is  the 
ionization  energy  of  a  hydrogen  atom;  oCH  *  0.666  X  10~24  cm3  is  the 
polarizability  of  a  hydrogen  atom;^(j  is  the  polarizability  of  the 
trap;  mQ  is  the  free  electron  mass;  and  m  is  the  effective  electron 
mass.  Lax^l  yggd  the  polarization  potential  well  to  model  capture 
of  electrons  and  holes  by  neutral  impurities.  Tasch  and  SatJ44^  used 
this  potential  well  to  model  the  observed  field  dependence  emission 
rate  of  Au  in  silicon.  Neutral  impurity  trap  such  as  AlGa°  has  the 
properties  of  a  polarization  potential  well.  In  the  presence  of  an 
electric  field,  the  total  potential  of  polarization  well  is  given  by 

VT(r)  *  -A<|/r4  -  qFr  cos(0)  (4.36) 
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4. 3. 4.1  Boole-Frenkel  effect 


The  minimum  potential  is  obtained  by  setting  dVT(r)/dr  =  0  at  r  = 
rfnax,  and  the  result  is  given  by: 


rmax  =  [4A1/cJFcos(W]1/5 
Eti  =  -1.649  A^IqFcostfl)}4/5 


(4.37) 


Note  that  the  emission  rate  is  field  dependent  for  0  <  0  <  ft/2,  and  is 
independent  of  electric  field  for  ft/2  <  0  Of.  Thus,  the  normalized  field 
enhanced  emission  rate  due  to  Poole- Frenkel  effect  can  be  written  as: 

eni/eno  =  exp(R)  (4.38) 

for  the  one  dimensional  case,  and 

en3/eno  "  0/2)  U+  <5/4)  f1  Rt1/4exp(R)dt)  (4.39) 

for  the  three  dimensional  case.  Here  R  =  1.649  A^5  (qF) 4//5 AT. 

4.3.4. 2  Phonon-assisted  tuoneliog  effect 

Substituting  Eq.(4.36)  into  Eq.(4.13)  and  letting  z  =  (Efcj  -  Eth)/kT 
yields  the  normalized  enhanced  emission  rate  due  to  the  phonon-assisted 
tunneling  effect: 

(E^- ^Eti  )/kT  *i/9  a 

ent/eno*f  exp{z  -[(8m  )  1/2/W  (-Aj/r4 

'o 

-  qFx+zkT)1,  2  dx}  dz  (4.40) 
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Numerical  integration  can  be  carried  out  by  assuming  x^  =  (A^  /  zkT) and 
x. 2  =  zkT/qF.  The  total  enhanced  emission  rate  is 

enhp  =  ®no  ^3  +  ®nt^  (4.41) 

Note  that  integrals  given  in  Eq.(4.39)  and  (4.40)  can  not  be  solved 
analytically.  The  numerical  calculations  of  the  enhanced  emission  rates  as 
a  function  of  the  electric  field  are  shown  in  Fig.4.5  for  Ec  -  0.35  eV  trap 
in  GaAs  at  300  K.  The  Poole- Frenkel  effect  is  dominant  when  electric  field 
is  lower  than  4  x  104  V/cm,  while  phonon-assisted  tunneling  effect  becomes 
important  for  electric  field  greater  than  4  x  104  V/cm. 

4.3.5  Dipole  potential  well 

The  dipole  potential  can  be  represented  by  a  trap  center  with  two 
oppositely  charged  ions,  each  with  charge  Zq. 

V(r)  =  r.  p/4/l*0*sr2  (4.42) 

where  p  is  the  dipole  moment,  p  =>  ZqRQ;  RQ  is  the  distance  between  two 
oppositely  charges.  Anti  site  pair  defect  (AsGa++GaAs  )  behaves  as  the 
dipole  potential  well  in  the  capture-emission  process.  As  is  shown  in 
Fig.4.6,  if  the  applied  electric  field  ,F,  forms  an  angle  0F  with  z  axis, 
and  is  polarized  along  the  z  direction,  then  the  total  potential  becomes: 

VT(r)  =-qPcos  @/47t€Qesr2-  qFr  si  n (  9p)  si  n  ©  cos  f 

-  qrF  cos($£0cos(e)  (4.43) 

4.3. 5.1  Poole-Frenkel  Effect 

The  mininum  potential  is  obtained  by  setting  dVT(r)/dr  =  0,  which 
yields: 


Normalization  Emissi 


'  Polarization  well 

2  Phonon  as  is  ted 

3  Poole  Frenkel 


P  (  104  V/cm) 


Flg*4*5  Nowalized  enhanced  electron  • 

lays 

sssr***  -“u«»  j? 


rmax  =  (p  cos  (g)/2  it  [F  sin^sinffl  cos($ 

+  Fcos(fijp)cos(^)  ]  J1/3  (4.44) 

,  and  the  change  in  barrier  height  is  given  by 

Eti(0,?>  =  -(3  x  2"2/3)  [pcos<0)/4 x e0er] 1/3 

x[Fsin(0p)sin(0)cos($  +  Fcos(6^)cos(0)]2/3  (4.45) 

Unlike  other  types  of  potential  wells,  the  dipole  potential  experiences  a 
barrier  lowering  which  is  field  dependent.  The  integral  falls  into  two 
cases  according  to  the  incident  angle  ,0p. 

1.  Casel:  0  <  <  k/2 

The  barrier  is  lowered  for  -  */2  <  §  <  Tt/2,  and  for  "X/2  <  $  <37C/2, 
0  <  9  <  tan”1  (-cot  (0F) /cos  (f)  ] . 

The  normalized  field  enhanced  emission  rate  is  thus  given  by: 

*/2  ?2 

en3/eno=  U/2*)  {  dj  {  sin(5)  exp[-/?  AEti  (  0  ,j>)  ]  d0 
' -iy 2  *  o 

37^2  tan-1 (-cotflp/cos (§) ) 

+  (1/270  f  df  (  s:n(  )exp[-£dEti  (0,4>)]d* 

1  7T/2  '  o 

3*/2 

+  (1/2X)  J^{l+[cot(^F)/cos(g)]2}"1/2  d#  (4.46) 

2.  Case  2:  7t/2  ^9^  <71 

The  barrier  is  increased  for  7t/2  <  £  <3  71/2,  and  for  -^2  <  $  <  n/2, 
0  <fl<  tan"1  [-cot (^p) /cos (5)1.  The  barrier  is  lowered  for  -^2  <  £<77/2 
,  and  for  tan"1[-cot(tfp)/cos($)]  <  0  <*/2. 


The  normalized  field  enhanced  emission  rate  for  this  case  is  given  by: 


en3/eno=l  -  (1/2*  )  /  U+  [cot(0p) /cos («#>)]  2 r1/2  d<t 

m  192 

+  (1/2  )/  dij'  sin<0)  exp[-<3  AEti«0,f)]d0  (4.47) 

-192  '  tan-1(-cot(5F)/cos($)) 

For  0F  *  0,  the©  dependence  disappears,  and  the  problem  can  be  solved 
analytically.  The  barrier  lowering  simplifies  to 

Eti  *  -1.9  (p /47T60€-r)1/3  F2/3  cos  Iff)  (4.48) 

,  and  the  normalized  Poole-Frenkel  effect  emission  rate  enhancement  is 
given  by 

en3/eno  =  d/R>  (exP(R)  -  U  (4.49) 

where  R  -  1.9  (p/4*^^) 1/3  F2/3AT. 

4.3. 5. 2  Ftnoao-nassisted  Tumeling  Effect 

For  the  potential  wells  studied  here,  the  one-dimensional  analysis  was 
done  for  0*0,  where  the  Poole-Frenkel  harrier  lowering  is  dominant.  For 
the  dipole  well,  the  maximum  barrier  lowering  occurs  at  9  =  <9m,  where  0  < 
6m  <9 F  and  ^  *  0.  Now,  differentiating  Eq.(4.46)  for  4Etj  with  respect 
to  6  at  ^*0,  and  setting  the  result  equal  to  zero,  we  obtain  a  solution 

foiGm 

$ra  *  tan-1  {(8  +  cos2(0p))1/2  -  3  cos(0p)  ]/2sin(0p) }  (4.50) 

Substituting  Eq.(4.50)  into  Eq.(4.44)  and  Eq.(4.44)  into  Eq.(4.13),  and 
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letting  z  =  (Etj  -  Efch)  /  kT,  the  field  enhanced  emission  rate  due  to 
phonon-ass i sted  tunneling  effect  can  be  expressed  as 


*nt/eno  =j 

'a 


(Etj-/lEti)AT 


exp{z  -  [{8m 


[-qPcos  {9m)/4neoesr* 


-qrF  sin(0F)sin(flj^-qrF  cos(dj)  cos  (0^+zkT]  ^%r}  dz  (4.51) 


Here  r^  =  (qPcosfi^^fQ^zkT) and  ^  =  zkT/qF[sin(^0  sin(0m)  +  cos(0p) 
cos(0m)].  The  normalized  field  enhanced  emission  rate  due  to  phonon- 
assisted  tunneling  effect  and  Poole-Frenkel  effect  can  be  calculated  by 
numerical  method.  The  results  are  shown  in  Fig.4.7.  Thus,  the  total  field 
enhanced  emission  rate  is 


enHD  *  ®no  (en3  +  *nt>  (4*52) 

Table  4.1  shows  the  plots  for  the  Coulombicpotential  well,  Dirac  well, 
square  well,  polarization  well  and  the  dipole  well.  Table  4.2  summarizes 
the  one-dimensional  and  three-dimensional  emission  rate  due  Poole-Frenkel 
effect  for  five  different  potential  wells.  Fig.4.8  shows  the  emission 
rate  vs  electric  field  for  the  Dirac  well,  square  well,  polarization  well 
and  the  Coulombic  well  as  compared  with  the  zero  electric  field  values. 
The  enhanced  emission  rates  are  nearly  identical  for  the  neutral  trap  with 
Dirac  well,  square  well  or  polarization  well.  However,  the  enhanced 
emission  rate  for  the  Coulombic  well  depends  very  strongly  on  the  electric 
field. 

4.4  Theoretical  Calculations  of  the  Nonezpooential  DLTS 
npopnnno  for  Different  Potential  Wells. 


From  the  analysis  of  field  enhanced  emission  rate,  the  electric  field 
dependent  DLTS  response  can  be  calculated.  This  is  discussed  as  follows: 


Normalization  Emission  Rate 


Normalization  Emission  Rate 


Consider  a  p+/n  abrupt  junction  diode,  the  electric  field  is  in  general 
spatial  dependent  within  the  depletion  region  and  can  be  expressed  by: 


F  =  W1  -  x/«) 


(4.53) 


where  x  is  the  distance  from  the  junction;  Fmax  is  the  maximum  electric 
field  occurred  at  the  metallurgical  junction  of  the  diode,  and  is  given  by 


Fmax  ^Ndw/  ^  o^r 

W  is  the  depletion  width  under  reverse  bias  condition. 
W  =  I2%€s  (Vbi  +  Vr)/qNd] 1/2 
=  LD[2(Vbi+Vr)AT  -  2] 1/2 
Lq  is  the  extrinsic  Debye  length  which  is  given  by 
Lp  *  [C^kT/tqftty]1/2 
vbi  ^  the  built-in  potential  which  reads: 


(4.54) 


(4.55) 


(4.56) 


vbi  =  (kT/ q)  In  (ND/n  j ) 


(4.57) 


n^  is  the  intrinsic  carrier  density, 
nj  =  (NcNv)1/2exp(-Eg/2kT) 

=4.9xl015(mdemdh/mo2)3/4T3//2  exp(-Eg/2kT)  (4.58) 

and  Eg  is  the  energy  bandgap,  which  is  a  function  of  temperature  ^  41]* 

Eg(T)  *  Eg(0)  T 2/(T  +r2)  (4.59) 

Eg(0)  *  1.519eV,  -  5.405  x  10_4  eV/K,  and  r2=  204  for  GaAs.  Electrons 
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which  are  located  in  the  region  between  WQ  (i.e.,  zero  bias  depletion  layer 
width)  and  W  (reverse  bias  depletion  layer  width)  will  be  emitted  into  the 
conduction  band  when  an  applied  bias  is  increased  from  0  to  -Vr,  as  is 
shown  in  Fig.4.9.  The  depletion  width  in  the  junction  space  charge  region 
is  first  divided  into  equally  spaced  small  segments.  In  each  small 
segment,  both  the  electron  trap  density  and  electric  field  are  assumed 
constant.  For  each  electric  field  (F^)  strength,  there  is  a  corresponding 
enhanced  emission  rate,  enj(Fj).  If  one  assumes  that  the  electron  emission 
is  exponential  transient  within  each  segment,  namely,  exp(enj t) ,  then  the 
total  emission  transients  in  the  depletion  layer  region  is  equal  to  the  sum 
of  the  individual  components,  which  can  be  expressed  by: 

<0 

Z  exp[-eni(Fj)t]  (4.60) 

Thus,  the  DLTS  signal  for  the  nonexponential  transient  can  be  expressed  by 

S(t)  expl-e^lF^t^  -  2  expt-eni  (Fjjtjl  (4.61) 

where  enj  is  the  electron  emission  rate  within  each  small  segment  of  the 
depletion  region,  as  is  shown  in  Eq.(4.14).  Values  of  en^  can  be 
calculated  for  different  potential  wells  using  the  equations  derived  above. 
From  Eq. (4.61),  it  is  noted  that  the  DLTS  spectral  response  is 
nonexponential.  Fig.4.10  illustrates  the  DLTS  response  for  the  EL2 
electron  trap  in  GaAs  calculated  from  Eq.(4.61)  for  different  potential 
wells  and  for  zero  electric  field.  Theoretical  calculations  showed  that  the 
location  of  DLTS  signal  peak  for  the  attractive  or  neutral  trap  is 
different  along  the  temperature  axis.  For  Coulombic  potential  well,  the 
DLTS  signal  peak  occurs  at  the  lowest  temperature,  while  for  the  neutral 
trap  with  Dirac  well,  square  well,  or  polarization  well  the  peak  of  the 
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DLTS  signal  occurs  nearly  at  the  same  temperature.  Fig.4.11  shows  the  DLTS 
spectral  response  for  the  EL2  electron  trap  for  the  cases  of  Coulombic  well 
with  single-  and  double-charge  state  as  well  as  the  experimental  results. 
Fig.4.12  illustates  the  DLTS  signal  for  the  EL2  electron  trap  in  GaAs  for 
a  Coulombic  well  with  a  double  charge  states  for  different  window  rates. 
By  comparison  of  the  theoretical  and  experimental  DLTS  spectral  response  vs 
temperature  will  enable  us  to  determine  the  potential  well  for  the  EL2 
electron  trap.  Our  results  showed  that  the  most  probable  potential  well  for 
the  EL2  trap  is  the  Coulombic  well  with  a  double  charge  states  (e.g., 
ASgg.n.  arsenic  anti  site  defect). 


1.  Coulombic  well 

2.  Dirac  well 


square  well 
polarization  well 


300 


350 
T  (K) 


400 


Fig. 4. 11  Calculated  DLTS  response  for  the  EL2-a  trap  in  GaAs,  assuming 
Coulanbic  well  with  (1)  double  charge  state  and  (2)  single¬ 
charge  state.  c_ 


i  Fig. 4. 12  Calculated  DLTS  response  for  the  EL2-a  trap  in  GaAs,  assuming 

*"  Coulombic  well  with  double  charge  state  for  five  different 

electron  anission  rates. 
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Table  4.1  Suimary  of  the  potential,  Poole-Frenkel  barrier 
for  five  different  potential  wells. 


V.  MODELLING  OP  NATIVE  POINT  DEFECTS  IN  GaAS 


Theoretical  modelling  of  native  point  defects  in  GaAs  is  presented  in 
this  chapter.  Based  on  chemical  reaction  and  principles  of  thermodynamics, 
expressions  for  the  equilibrium  defect  concentration  as  functions  of 
temperature  and  arsenic  pressure  during  crystal  growth  are  derived. 
Thermal  kinetic  equations  are  then  employed  to  predict  the  possible  native 
defects  in  GaAs  after  crystal  cooling.  For  GaAs  grown  under  As-rich  or 
high  arsenic  pressure  condition,  it  is  shown  that  several  native  defects 
such  as  gallium  vacancy  (V^),  arsenic  interstitial  (ASj),  arsenic  anti  site 
(Asgg),  arsenic  anti  site-plus-arsenic  vacancy  (V^^ASqq)  and  their  complexes 
may  be  observed  in  GaAs  material. 


5.1  Introduction 

GaAs  specimens  grown  by  various  techniques  such  as  Liquid 
Encapsulation  Czockraski  (LEC) ,  Vapor  Phase  Epitaxy  (VPE),  Liquid  Phase 
Epitaxy  (LPE),  Organometalllic  Chemical  Vapor  Deposition  (MOCVD),  and 
Molecular  Beam  Epitaxy  (MBE)  are  known  likely  to  produce  different  defect 
properties  in  terms  of  energy  levels  within  the  bandgap  of  a  semiconductor 
due  to  differences  in  the  native  defects  and  trace  impurities  or  impurity 
complexes.  In  each  of  these  techniques,  the  growth  temperature,  growth 
pressure,  growth  phase,  and  cooling  rate  are  usually  different.  For 
example,  GaAs  grown  by  LPE  technique  from  a  gallium  melt  is  expected  to  be 
low  in  gallium  vacancy  defects  and  high  in  arsenic  vacancy.  The  LPE  GaAs 
usually  contains  hole  traps  with  energy  level  of  Ev  +  0.71  eV  (B-center), 
while  the  VPE  GaAs  grown  in  high  arsenic  pressure  condition  always  contains 
an  electron  trap  with  energy  of  Ec  -  0.83  eV  (EL2  level). 
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In  this  chapter,  the  chemical- thermodynamic  principles  and  thermal 
kinetic  equations  are  introduced  to  model  the  grown-in  defects  in  GaAs. 
Using  chemical-thermodynamic  principles  to  analyze  point  defects  in  GaAs 
was  first  proposed  by  Logan  and  Hurle  1^2]  m  They  considered  the  shallow 
vacancy  as  well  as  shallow  interstitial  levels.  Bublik  modified  their 
model,  and  applied  it  to  the  deep  vacancy  level.  More  recently,  Hurle 
used  Frenkel  defect  instead  of  Schottky  defect  for  the  arsenic  vacancy  and 
arsenic  interstitial  defects.  During  crystal  growth,  both  VGa  and  VAs 
Schottky  vacancy  pairs  are  important  defects  for  the  undoped  GaAs.  In  the 
present  work,  we  use  vGa  and  VAS  Schottky  vacancy  pairs  with  As-Frenkel 
defect  to  calculate  the  vacancy  and  interstitial  concentrations  during 
crystal  growth.  During  crystal  cooling  from  high  temperature  to  room 
temperature,  thermal  kinetics  is  expected  to  be  the  dominant  mechanism. 
Since  the  formation  energy  for  anti  site  defect  is  very  low  in  the  after¬ 
growth  condition  (VGa  +  Asj  ■  AsGa,  H  *  0.35  eV) ,  the  anti  site  defect  is 
likely  to  be  the  dominant  defect  after  crystal  growth,  as  was  proposed  by 
Van  Vechten  [®*1. 

From  the  theoretical  analysis  of  GaAs  material  grown  under  As-rich  or 
high  arsenic  pressure  condition  (e.g.,  LEC,  VPE,  and  MOCVD  techniques), 
defects  such  as  arsenic  interstitial  (As^),  gallium  vacancy  (Vgg),  arsenic 
anti  site  (As^) ,  arsenic  anti  site-plus-arsenic  vacancy  (V^gAs^),  and  their 
complexes  are  the  possible  defects  in  undoped  GaAs,  while,  gallium  vacancy 
(V^),  arsenic  vacancy  (VAs)  and  their  complexes  are  the  possible  grown  -in 
defects  in  both  the  MBE  and  LPE  GaAs  under  low  arsenic  pressure  condition. 

Section  5.2  deals  with  the  theoretical  calculations  of  vacancy  and 
interstitial  defects  in  GaAs  during  crystal  growth.  Thermal  kinetic 
equations  are  described  in  section  5.3.  Section  5.4  discusses  the  possible 


grown-jn  point  defects  in  GaAs  for  the  As-rich  or  high  arsenic  pressure 
case.  The  possible  grown-in  point  defects  in  GaAs  for  the  Ga-rich  or  low 
arsenic  pressure  cases  are  depicted  in  section  5.5. 

5.2  Theoretical  Calculations  of  Vacancy  and  Interstitial  in  Ohdnped  GaAs 

In  this  section,  the  chemical-thermodynamic  principles  162—69]  ^ 
employed  to  derive  expressions  for  thermal  equilibrium  defect  concentration 
as  functions  of  temperature  and  arsenic  pressure  in  GaAs.  Hurle  [63] 
calculated  the  point  defects  by  considering  only  As-vacancy  and  As- 
interstitial  defects  in  undoped  GaAs.  However,  it  is  known  that  Ga- vacancy 
related  defects  may  play  an  important  role  in  anti  site  formation. 

In  this  chapter,  we  consider  several  types  of  point  defects  such  as 
arsenic  monovacancy,  VAg,  positively  charged  arsenic  monovacancy,  VAg+, 
arsenic  interstitial,  As^,  positively  charged  arsenic  interstitial,  ASj+, 
gallium  monovacancy,  v^,  and  negatively  charged  gallium  monovacancy  V,^-. 

To  deal  with  the  problem  of  native  point  defects  in  GaAs,  one  can 
write  down  the  kinetic  reaction  equations  for  formation  of  each  type  of 
defects  as  well  as  formation  of  electrons  and  holes  in  the  crystal.  There 
is  an  additional  reaction  equation  which  represents  the  transfer  of  atoms 
between  gas  and  solid  phases.  For  each  of  these  reactions,  there  is  a  mass 
action  which  applies  in  equilibrium.  The  mass  action  can  be  written  in 
terms  of  concentration.  To  these  mass  action  equations  one  adds  the 
condition  of  charge-neutrality,  and  the  resulting  set  of  equations  can  then 
be  solved  as  a  function  of  arsenic  partial  pressure.  Considering  the 
defects  cited  above,  one  has  the  following  reaction  equations. 

*  Vi  *  *>i  +  <5-» 

Vfta  *  W  *  e  '  <5-» 
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Asj  =  Asj  +  e 


0  =  e“  +  h 


(5.4) 


(l/2)As2(g)  +  Wi  =  ASj 


(5.5) 


0  "  vGa  +  vAs 


(5.6) 


VGa  «  vGa"  +  h 


(5.7) 


Eqs.(5.1)  through  (5.4)  represent  the  reactions  for  forming  the  ionized 
arsenic  Frenkel  defects.  Eq.(5.5)  denotes  the  transfer  of  arsenic  atoms 
between  solid  and  gas  phases.  Eq.(5.6)  shows  the  formation  of  Schottky 
pairs.  Eq.(5.7)  represents  the  ionization  of  a  Ga- vacancy.  The  mass 
action  relationships  corresponding  to  the  above  reactions  are  given  as 
follows: 162,63] 


Kfa  "  [ASi]  [Vas] 


=  2.92  x  106  exp (-4. 845 AT) 


(5.8) 


=  n  [V+]  /  [V. J 


«  442  exp (-0.27 AT) 


Kai  *  n  [Vj+]  /  (Asj] 


(5.9) 


*  4.9  x  10^  x  T^2  exp (-0.4 AT) 


(5.10) 


*  nP 


*  1  x  10-12  T2  exp  (-1.62 AT) 
KM2i  ‘  »»jl  PAS2'1/2 


(5.11) 


*  16.4  exp (-1.125 AT) 


(5.12) 
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«s  =  fvGal  tVAs^ 


=  3.286xl04  exp (-3. 6 AT)  (5.13) 

V  "  p  tvGa‘l  /  tV ca] 

=  3.7  x  10-8  T3/2  exp (-0.66 AT)  (5.14) 

Square  brackets  in  the  above  equations  denote  the  concentration.  The  As^ 
represents  arsenic  atom  at  the  arsenic  site,  and  is  taken  as  unity.  The 
partial  pressure  of  AS2  in  the  gas  phase  is  denoted  by  PAs2.  The 
equilibrium  constants  appearing  on  the  right  handside  of  Eqs.(5.8)  through 
(5.14)  have  the  general  form: 

K  *  exp(  SA)  exp(-  H/kT)  =  KQ  exp(-  HAT)  (5.15) 

where  S  and  H  are  the  entropy  and  enthalpy  changes  for  each  reaction, 
respectively.  The  charge  neutrality  condition  is  obtained  from  Poisson 
equation: 

n  +  tV^H  *  p  +  [Vas+]  +  [ASj+]  (5.16) 

Solving  Eqs.(5.8)  through  (5.16)  one  obtains  an  expression  for  n2  as: 

"2  "  +  <KavKfa/KAe2i>  pAs2_1/2  +  KaiKss2ipAs21/2) 

/<1+VKsKAs2ipAs21/2/KfaKcv'  <5-17> 

Thus,  the  electron  concentration,  n,  can  be  calculated  from  Eq.(5.17),  and 
other  defect  density  can  also  be  deduced  from  n  via  Eqs.(5.8)  to  (5.16). 
Fig.5.1  shows  the  defect  concentration  vs  PAs2  for  T  =  1000  K.  Note  that 
Ga  vacancy  is  the  dominant  point  defect  in  the  entire  As2  pressure  range 
shown.  The  LPE  Ga  As  is  usually  grown  in  the  lower  As2  pressure  range  with 
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are  the  dominant 


pAs2  *  l0-8  atm*  and  defects  such  as  VGa,  VGa~,  and  VAs+ 
defects.  For  VPE  GaAs  with  a  corresponding  partial  pressure  of  1 0“4  atm, 
defects  such  as  VGa,  VGa“,  VAs+,  ASj+,  and  ASj  are  the  dominant  defects. 
In  the  pressure  range  between  LPE  and  VPE  growth,  the  concentration  of  VGa" 
is  equal  to  the  concentration  of  VAS+*  Figure  5.1  also  shows  that  the 
grown- in  defects  for  the  VPE  and  MOCVD  techniques  are  much  more  than  LPE 
and  MBE  growth  techniques.  Fig.5.2  shows  defect  density  as  a  function  of 
growth  temperature  for  PAs2  =  5  x  10  atm.  The  concentration  of  VGa 
monovacancy  and  VAS+  monovacancy  increases  with  increasing  growth 
temperature,  but  the  density  of  ASj  +  mono-interstitial  decreases  with 
increasing  growth  temperature. 

5.3  Tbeonl  Kinetics  after  Crystal  Growth: 

In  general,  defects  will  migrate  as  the  crystal  cools  down  from  the 
growth  temperature.  Thus,  one  would  expect  defect  concentrations  to  reach  a 
new  equilibrium  condition  at  lower  temperature  The  enthalpy  of  a 

single  vacancy  migration  in  GaAs  is  1.6  eV  [64, 72]^  For  examples,  if  the 
jump  attempt  frequency  is  equal  to  Debye  frequency,  then  the  jump  rate  will 
be  4  x  1 07  s_1  at  1400  K  and  5  x  104  s_1  at  1000  K.  For  example,  for  a 
typical  crystal  to  cool  down  from  1050  K  to  950  K  in  10  min.,  simple 
vacancy  would  pass  through  more  than  10®  lattice  sites  in  that  period.  The 
migration  of  vacancy  will  produce  anti  site  defects,  anti  site  pairs  or 
impurity  complexes.  In  LPE  grown  GaAs  samples,  we  have  observed  that 
samples  with  faster  cooling  rate  {1°  C/min.)  would  produce  more  defects 
than  slower  cooling  rate  (0.4°  C/min.).  Most  of  the  simple  vacancies  are 
ionized,  and  may  encounter  with  other  defects  to  form  bound  complexes.  At 
room  temperature,  almost  all  the  vacancies  present  are  likely  to  interact 
with  other  defects  to  form  complexes.  To  consider  the  ultimate  fate  of  a 
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simple  point  defect  introduced  during  crystal  growth,  several  types  of 
defects  and  defect  complexes  should  be  considered.  It  seems  likely  that 
the  most  common  defect  complexes  will  be  those  which  have  no  net  charge 
because  these  would  have  the  greater  binding  energy.  For  example,  a 
negatively  charged  gallium  vacancy  and  a  positively  charged  arsenic  vacancy 
may  produce  a  neutral  gallium-arsenic  di vacancy.  This  can  be  expressed  by: 


V  +  V  *  vGa'vte 


(5.18) 


In  addition,  the  anti  si  te-di  vacancy  complexes  may  be  formed  by  the 
reaction  equations  shown  below: 


and 


2Vfts  +  Ga 


-  _  tr  +, 


AS 


*><*  *  *=Ga  ‘  vGa'AsGa  ^  vGa' 


(5.19) 

(5.20) 


The  anti  site  pair  defect  complexes  can  be  written  as: 


G®as  +  ^Ga  *  ^As  ^ 


Ga 


(5.21) 


The  interaction  of  Vfts+  (VGa“)  with  a  single  acceptor  (or  donor) 
impurity.  A”  (D+),  can  be  written  as: 


A  +VAs  3AAs 


Dt  +  V, 


Ga  uGa 


(5.22) 

(5.23) 


which  should  have  about  the  same  binding  energy  as  the  gallium-arsenic 
divacancy  complex  given  by  Eq.(5.18).  From  Fig. 5.1,  it  is  noted  that  the 
concentration  of  VGa  is  very  high,  and  the  following  reactions  are 
prevailed: 


(5.24) 


'As 


+  Vca  +  e"  =  VnaV 


Ga  As 


VQa  +  D+  "  D< 


'Ga 


(5.25) 


In  general,  the  concentration  of  positively  charged  and  negatively 
charged  simple  point  defects  introduced  during  crystal  growth  are  not 
equal.  Therefore,  the  concentration  of  one  type  of  simple  defects  may  be 
exhausted  as  shown  by  the  reaction  equations  (5.18)  to  (5.25).  This  can  be 
illustrated  by  the  following  kinetic  equations: 

(1)  Arsenic  antisite  can  be  formed  by  ASj+,  ASj,  V^-,  and  vGa‘ 


ASj+  +  =  AsGa++  +  2e“ 

(5.26) 

Asi+  +  vGa  “  ^Ga4*  +  e' 

(5.27) 

^i  +  VGa~  ■  +  3e" 

(5.28) 

ASj  +  =  ASgg^  +  2e“ 

(5.29) 

(2) 

Antisite-vacancy  complexes  can  be  expressed  by: 

vGa~  ♦  + 

(5.30) 

vGa  +  *  “oa VAS~+  * 

(5.31) 

vAs+  +  C^Ga  *  *  S>AsvGa_ 

(5.32) 

which  will  occur  everytime  a  vacancy  migrates  to  a  nearest  neighbor 

(3) 

Impur j ty- vacancy  complexes  can  be  written  as: 

and 

vGa  +  AAs  *  AAsvGa 

(5.33) 

vGa  +  °Ga+  +  e  3  °GavGa 

(5.34) 
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Eq. (5.33)  and  (5.34)  are  the  interaction  of  VGa~  with  a  acceptor  (donor) 
substitutional  site.  For  examples,  Te^V^"  defect  complex  may  be  formed  in 
Te-doped  GaAs,  while,  Sn^Vog-  and  GeQgVgg"  defect  complexes  may  be  formed 
in  Sn-doped  and  Ge-doped  GaAs,  respectively. 

The  VAg+  may  migrate  to  a  donor  substitutional  site,  and  the  reaction 
equation  is  given  by: 

W  *  Cfe.  ■  °eavAs*  <5-35> 

Obviously,  such  reactions  act  as  to  compensate  the  dopant  and  to  prevent 
the  Fermi -level  from  approaching  either  band  edge  during  crystal  cooling. 

5-4  Hie  Possible  Growo-in  Defects  in  GaAs  Older  As-rich 
or  High  Arsenic  Pressure  Condition 

Applying  defect  model  to  the  high  arsenic  pressure  conditions,  the 
VGa'  vGa"'  vAs+'  ASj+,  and  As j  are  found  to  be  the  dominant  defects  when 
the  arsine  partial  pressure,  PAg2,  is  greater  than  10-5  atm  (e.g.,  GaAs 
grown  by  LEC,  VPE,  and  MOCVD  techniques);  this  is  shown  in  Fig.5.1.  The 
migration  of  VGa'  vGa~'  and  vAs  defects  may  result  in  forming  complexes 
given  by  Eqs.(5.18)  to  (5.20)  and  Eqs.{5.22)  to  (5.35)  after  crystal 
growth.  If  we  neglect  the  neutral  defect  complexes,  then  native  defects 
such  as  AsGaVAs+++,  AsGa++,  VGa",  GaAgVQa  ,  VAg+,  ASj+,  VGa,  ASj ,  and 
impurity  complexes  Dog"*-,  AAgVGa~,  DQgV^”,  DGaVAg+  are  the  possible  defects 
which  may  be  performed  in  high  arsenic  pressure  conditions.  In  As-rich 
condition,  the  concentration  of  VAg+  may  be  decreased,  Asj  +  may  be 
increased,  and  AsGaVAg+++,  AsGa++,  VGa”,  ASj+,  VGa,  Asj,  and  impurity 
complexes  are  the  possible  defects. 


S.S  The  Possible  Grown- in  point  Defects  in  GaAs  far  the 
Ge-rich  or  Low  Arsenic  Pressure  Case. 

Apply  defect  model  to  the  low  arsenic  pressure  condition,  VGa'  vGa”' 

and  VAs+  are  found  to  be  the  dominant  point  defects  when  PAs2  is  less  than 

10“  7  atm;  this  is  the  case  for  the  LPE  and  MBE  grown  GaAs.  The  vacancy 

migration  after  crystal  growth  will  result  in  the  formation  of  complexes 

given  by  Eqs.(5.18)  to  (5.20),  Eqs.(5.22)  to  (5.24),  and  Eqs.(5.30)  to 

(5.35).  If  we  neglect  the  neutral  defect  complexes,  then  native  defects 

such  as  vGa"'  vas+'  <*a.vgT"  ”,  AsGa^As+++,  Vgg  and  impurity  complexes  such 

as  DGa+,  AASVGa_/  DGaVGa~'  °GaVAs+  are  the  P°ssible  defects  in  low  arsenic 

pressure  conditions.  The  concentration  of  vGa"  is  almost  equal  to  that  of 

VAs+'  as  is  shown  in  Fig.5.1,  and  the  binding  energy  for  the  gallium- 

arsenic  divacancy  is  very  high,  as  given  in  Eq.(5.18).  Thus,  high  purity 

GaAs  can  be  grown  in  this  case  under  optimum  cooling  condition.  In  Ga-rich 

and  low  arsenic  pressure  conditions,  the  concentration  of  and  may 

be  decreased,  gallium  interstitial  may  be  reacted  with  VAs+  to  form  GaAg  . 

Thus  the  possible  defects  are  VAg+,  Ga^  ,  and  Ga^V^  in  undoped  GaAs 

grown  under  Ga-rich  condition. 

5.6  Suaaury  and  Conclusions 

In  this  chapter,  we  presented  a  new  defect  model  for  predicting  native 
point  defects  in  GaAs.  During  crystal  growth,  the  chemi cal- thermodynamic 
principles  are  used  to  derive  the  density  of  vacancy  and  interstitial 
defects  under  thermal  equilibrium  condition.  After  crystal  growth,  thermal 
kinetic  equations  are  employed  to  predict  the  anti  site  defects  and  impurity 
complexes.  Conclusions  are  listed  as  follows: 

(1)  High  purity  GaAs  material  can  be  grown  for  the  low  arsenic  pressure 
case  under  optimum  cooling  condition. 


(2)  GaAs  grown  under  high  arsenic  pressure  condition  will  produce 

more  native  point  defects  than  under  low  arsenic  pressure  condition. 

(3)  Native  defects  such  as  AsGa,  AsGaVAg,  VGa,  and  Asd  and  impurity 

complexes  such  as  D ,  AAsVGa~,  D^V^-  are  the  possible  defects  for 
GaAs  layers  grown  under  the  As-rich  or  high  arsenic  pressure  condition 

(4)  Native  defects  such  as  VGa,  VAg,  Ga AsVGa ,  AsGaVAs  and  impurity 
complexes  such  as  DGa,  AAgVGa,  DGaVGa,  DGaVAg  are  the  possible 
defects  for  GaAs  grown  under  low  arsenic  pressure  condition. 

(5)  Native  defects  such  as  VAg,  GaAg,  GaAgVGa  and  impurity  complexes 
such  as  DGaVAs+  are  the  possible  defects  for  GaAs  grown  under  Ga-rich 
and  low  arsenic  pressure  conditions. 

(6)  Arsenic  anti  site  (AsGa)  defect  can  only  be  observed  in  GaAs  grown 
under  As-rich  or  high  arsenic  pressure  conditions;  this  defect  can  not 
be  produced  under  low  arsenic  pressure  and  Ga-rich  conditions. 


VI.  IK  PHYSICAL  ORIGINS  OF  EL2  QJCnON  TRAP  IN  GaAs. 

(THEORETICAL  AND  EXPERIMBITAL  EVIDOCE) 

The  activation  energy  for  the  EL2  electron  trap  in  GaAs  reported  in 
the  literature  ranges  from  Ec-0.76eV  to  Ec-0.83eV.  The  physical  origin  of 
this  trap  is  a  subject  of  great  interests  in  recent  years.  A  large  number 
of  papers  has  been  devoted  to  this  subject.  Based  on  our  theoretical  model 
and  experimental  results,  we  found  that  EL2  level  is  formed  by  two  electron 
traps.  One  is  identified  as  the  EL 2a  (Ec-0.83eV)  electron  trap,  the  other 
is  denoted  as  the  EL2b  (Ec-0.76eV)  trap.  The  physical  origin  of  EL2a  level 
is  attributed  to  the  arsenic  anti  site  (As^)  defect,  whereas,  the  physical 
origin  of  the  EL2b  level  is  due  to  arsenic  anti  site-plus-arsenic  vacancy 
(AsGavAs)  <^e^ect  complex. 

Based  on  the  native  defect  model  presented  in  chapter  5,  modelling  of 
the  EL2  electron  trap  is  derived  in  GaAs  material  grown  under  As-rich  or 
high  arsenic  pressure  condition.  Furthermore,  from  calculations  of  the 
trap  density  for  the  MOCVD  and  VPE  grown  GaAs  with  As  to  Ga  mole  fraction 
ratio  greater  than  one,  it  is  found  that  the  density  of  EL2a  trap  is 
proportional  to  the  mole  fraction  ratio  of  (r-1)^^;  whereas,  the  density 
of  EL2b  trap  is  proportional  to  (r-1)1/4,  where  r  *  [As]/[Ga], 

Theoretical  calculations  of  the  nonexponential  DLTS  signal  for  the 
EL2a  trap  revealed  that  the  best  fitted  potential  well  for  this  charged 
center  was  due  to  a  Coulombic  potential  well  with  double  charge  state  (ie., 
AsGa++).  This  was  verified  by  taking  into  account  the  Poole- Frenkel  effect 
and  phonon-assisted  tunneling  effect  of  the  electric  field  enhanced 
electron  emission  rates  data  for  this  trap. 
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To  pcoviode  the  experimental  evidence  for  the  proposed  physical 
origins  of  the  EL2  electron  traps  desi cussed  above,  DLTS  measurements  were 
performed  on  VPE  GaAs  layers  grown  on  different  orientations  (e.g.,  (100), 
(211A),  (211B))  and  LEC  GaAs  layers  annealed  in  hydrogen  gas  at  different 
temperatures  (e.g.,  200,  300,  500°C).  From  the  results  of  our  annealing 
study  of  the  LEC  grown  GaAs,  it  was  found  that  EL2b  trap  would  disappear, 
and  EL2a  trap  would  emerge  in  the  DLTS  scan  as  the  annealing  temperature 
increases  to  500°C.  This  result  may  be  interpreted  by  the  model  of  EL2 
electron  trap  to  be  presented  in  this  chapter.  The  experimental  evidence 
which  may  be  used  to  support  the  modelling  of  EL2  will  be  presented  in 
chapter  VII. 

Section  6.1  reviews  the  EL2  electron  trap  in  GaAs.  Theoretical  model 
of  EL2  electron  trap  is  described  in  section  6.2.  Section  6.3  discusses 
the  method  of  determining  the  potential  wells  for  the  EL2  trap  levels  from 
analyzing  the  field  enhanced  emission  rates  deduced  from  the  nonexponential 
DLTS  data.  Summary  and  conclusions  are  given  in  section  6.4. 

6.1.  Review  of  the  02  Electron  Trap  in  GaAs. 

Activation  energy  of  EL2  electron  trap  covers  the  energy  range  from 
Ec-0.75eV  to  Ec-0.83eV,  as  reported  by  many  previous  investigators^73-91^. 
This  trap  has  been  observed  in  GaAs  grown  by  the  Bridgmann^77'7^ ,  LEC^79- 
85^,  VPE^86'87^,  and  MOCVD t88-91)  methods,  as  well  as  high  temperature 
heat-treated  GaAs  samples f9^l.  However,  this  level  was  not  observed  in  the 
LPE^93^  and  MBE^94^  grown  GaAs  epitaxial  materials.  Recent  studies  of  the 
LEC  bulk  grown  GaAs  reported  by  Taniguchi  et  al  t74^  have  found  that  Ec- 
0.77eV  electron  trap  exists  in  the  front  section  of  LEC  GaAs  ingot,  whereas 
,  Ec-0.82eV  level  is  the  dominant  trap  level  appeared  in  tail  section  of 
the  i ngot.Thesame  result  was  reported  in  MOCVD  grown  GaAs  epilayersby 


Watanabe  et  al  ^5]^  They  found  that  EL2a  (Ec-0.83eV)  was  the  dominant 
trap  level  for  GaAs  grown  at  720  to  740°C,  while  EL2b  (Ec-0.76eV)  was  the 
dominant  trap  level  for  GaAs  grown  at  630  to  660°C.  In  their  annealing 
studies.  Day  et  al  ^84^  reported  that  EL2b  was  observed  in  as-grown  MBE  n- 
GaAs  material;  the  results  further  revealed  that  this  level  can  be  annealed 
out  at  800°C  for  1/2  hour  or  at  700°C  for  1  hour.  Whereas,  EL2a  was  not 
observed  in  the  as-grown  MBE  n-GaAs.  The  EL2a  level  can  be  created  by 
thermal  annealing  process  and  its  trap  concentration  can  be  enhanced  by 
high  temperature  annealing.  These  results  indicate  that  EL2  elctron  trap 
may  be  due  to  two  different  trap  levels.  One  of  them  is  the  EL2a  (Ec- 
0.83eV)  electron  trap  and  the  other  is  0L2b  (Ec-0.76eV)  electron  trap.  Our 
model  and  experimental  evidence  have  shown  that  these  two  electron  traps 
have  different  physical  origins. 

Studies  of  the  LEC  grown  GaAs  (79'95]  ^ave  shown  that  n-type,  (S.I.) 
material  can  be  grown  only  from  melts  above  a  critical  As  composition,  and 
the  EL2a  level  was  native  defect  observed  in  this  material.  Ga-rich  melts 
were  found  to  yield  p-type,  low  resistivity  GaAs  crystal.  Ta  et  al^7^ 
reported  that  EL2a  level  was  observed  in  As-rich  GaAs  material.  GaAs 
samples  prepared  by  LPE  method  are  grown  from  a  Ga-rich  solution,  whereas, 
WE  GaAs  samples  are  commonly  prepared  in  an  As-rich  gas  ambient,  in  which 
growth  rate  and  surface  morphology  are  optimized^].  The  EL2a  level  was 
not  observed  in  as-grown  MBE  n-GaAs  material  ^94^.  However,  it  was  found 
that  EL2a  level  can  be  created  and  its  trap  density  can  be  enhanced  by  high 
temperature  annealing  (above  500°  C)^4^.  There  is  clear  evidence  that  the 
EL2a  trap  can  only  be  observed  in  GaAs  grown  under  As-rich  or  high  arsenic 
pressure  case.  The  high  arsenic  pressure  will  enhance  the  formation  of 
vGa'  and  As-rich  case  will  increase  the  number  of  As^89^.  Thus,  the  EL2a 
level  is  associated  with  gallium  vacancy,  arsenic  interstitial  or  their 


complexes  ,  and  is  not  related  to  oxygen  impurity  complex  ^ 97 

In  undoped  GaAs  material,  Watanabe  et  al^5,89]  have  shown  that  the 
concentration  of  EL2b  level  is  proportional  to  ([AsHjJ/fTMG])1/4.  However, 
the  concentration  of  EL2a  level  is  proportional  to  ( [  AsH-j] / [TMG] ) They 
also  found  that  densities  of  SL2a  and  EL2b  level  would  increase  with 
increasing  growth  temperature  for  the  MOCVD  grown  GaAs.  Bhattacharya  et 
ai[87l  observed  a  linear  dependence  of  the  EL 2a  trap  density  on  [As]/[Ga] 
ratio  in  MOCVD  GaAs.  Lagowski  et  al^78^  have  found  that  the  concentration 
of  EL2a  level  was  increased  with  increasing  As  pressure  during  Bridgemann 
bulk  growth,  while.  Miller  et  al^86^  found  that  the  density  of  EL2a  level 
was  increased  with  increasing  [As^l/tGaCl]  ratio  in  VPE  grown  GaAs.  Li  et 
al  also  found  that  the  density  of  EL2a  level  was  decreased  with 
decreasing  [AsCl3]/[Ga]  ratio  in  the  Ga-rich  VPE  GaAs. 

In  the  S-doped  GaAs,  Watanabe  et  al^75'89^  found  that  density  of  EL2b 

level  decreases  with  increasing  dopant  concentration  of  surfur  impurity. 

The  density  of  EL2a  level  was  found  to  decrease  with  increasing 

concentration  of  shallow  donor  dopants  {Si,  Se,  Te)  as  was  observed  by 

Lagowski  et  al t78^ .  Donor  concentration  above  a  threshold  value  (1  x  10^7 

cm--*)  led  to  a  rapid  elimination  of  EL2  trap^78^.  This  is  consistent  with 

our  observation ^ 98 ^  in  which  no  EL2  trap  level  was  detected  in  MOCVD  GaAs 
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with  Sn  dopant  density  higher  than  3  x  10  cm  .  If  group  VI  (Se,  Te,  S) 

elements  occupy  the  As  lattice  site  and  combined  with  gallium  vacancy  to 

form  anti  si  te-vacancy  (A^V^-)  complexes,  then  the  concentration  of  VGa 

will  be  decreased.  Thus,  the  probability  of  forming  an  arsenic  antisite 

(As defect  will  be  greatly  reduced  t87l.  Johnson  et  al  used 

photo 1 urn inescene  upconversion  method  to  observe  the  EL2a  level,  and 

concluded  that  this  level  is  due  to  arsenic  anti  site  defect.  The  fact  that 

EL2a  level  is  a  donor  type  defect  was  also  supported  by  the  observed  field 

dependence  of  emission  rates^48'7^.  This  is  consistent  with  the  modelling 

70 


of  native  point  defect  described  in  the  previous  chapter.  We  conclude 
that  arsenic  anti  site  defect  is  the  only  grown-in  defect  which  is  observed 
in  As-rich  or  high  arsenic  pressure  case,  but  is  not  observed  in  Ga-rich 
and  low  arsenic  pressure  case.  From  the  experimental  evidence,  it  can  be 
shown  that  EL2a  trap  is  due  to  anti  site  defect,  AsGa++,  formed  during  the 
post-grown  cooling  154,99] ,  as  ^  discussed  further  in  next  section. 

6.2.  Theoretical  Modelling  of  the  EL2  Electron  Trap  in  GaAs. 

We  shall  next  present  a  new  model  for  explaining  the  EL2  level  vs 
different  [As]/[Ga]  ratio  for  the  MOCVD  and  VPE  grown  GaAs  samples. 
Assuming  that  the  mole  fraction  ratio  of  [As]  to  [Ga]  is  equal  to  r,  it  can 
be  shown  that  for  r  >  1,  the  concentration  of  EL2a  level  is  proportional  to 
(r-1)1/2,  and  the  concentration  of  EL2b  level  is  proportional  to  (r-1)1/4. 
In  this  section,  the  kinetics  of  EL2  formation  in  GaAs  are  described.  In 
the  growth  process,  vacancies,  interstitials  and  anti  sites  are  formed  in 
high  temperature  thermal  equilibrium.  In  the  cooling  process,  vacancies 
migrate  to  form  EL2a  and  EL2b  electron  trap.  In  the  annealing  process, 
EL2b  trap  level  is  transferred  to  EL2a  trap  level  after  thermal  annealing. 
These  three  different  processes  of  defect  formation  in  GaAs  are  discussed 
next. 

6.2.1.  Growth  process: 

Defect  formation  under  thermal  equilibrium  for  the  case  of  r  >  1 
(a)  MOCVD  grown  GaAs  epitaxial  material: 

Arsine  (AsH3)  and  trimethylgallium  (TMG)  were  used  as  sources  for  As 
and  Ga  in  MOCVD  grown  GaAs.  Let  the  mole  fraction  ratio  of  arsine  to  TMG 
be  r,  then: 

[AsH3]/[Ga(GH3)3]  *  r  (6.1) 
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If  we  assume  [Ga  (013)3]  equal  to  1,  then  [ASH3]  will  he  equal  to  r, 
where  the  square  bracket  in  Eq.(6.1)  represents  the  mole  fraction  of  arsine 
and  TMG  gas.  In  the  crystal  growth  arsine  will  decompose  into  As4  (g)  or 
As2 (g) ,  depending  on  the  growth  temperature.  In  general,  arsine  will 
decompose  into  As4(g)  in  the  epitaxial  growth  (growth  temperature  below 
1000  °c),  while  arsine  will  decompose  into  As2  (g)  in  the  melt  growth  in 
which  the  growth  temperature  is  above  1400  °C^64^.  The  growth  temperature 
for  MOCVD  process  is  usually  below  1000  °C,  and  thus  the  reaction  of  arsine 
can  be  written  as: 

AsH3(g)  =  (l/4)As4(g)  +  3/2  H2(g)  (6.2) 

The  chan i cal  reaction  of  As4(g)  and  TMG  can  be  described  by: 

Ga (CH3) 3  +  (r/4)As4(g)  +  (3/2)H2(g) 

=  GaAs(s)  +  [(r-l)/4]As4(g)  +3CH4(g)  (6.3) 

If  GaAs(g)  is  completely  deposited  on  surface  of  GaAs  substrate,  then  the 
extra  As4(g)  will  be  decomposed  into  AsAs(s)  or  As^ ,  by  the  following 


reactions^®®! : 

(1/4)  As4  (g) 

=  ASas<s) 

+  VGa 

(6.4) 

(l/4)As4(g) 

*  Asi 

(6.5) 

(1/4) As4 (g) 

*  ^as^ 

+  ^Ga 

(6.6) 

Eqs.(6.4)  and  (6.5)  may  be  obtained  from  Eqs.(5.1)  through  Eq.(5.7).  These 
two  equations  are  the  dominant  reactions  during  crystal  growth.  The 
reaction  of  Eq.(6.6)  can  be  neglected  compared  to  Eqs.(6.4)  and  (6.5).  The 
law  of  mass  action  for  Eqs.(6.4)  through  (6.6)  is  given  by: 
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Fig. 6.1  As^g  or  As  atoms  are  deposited  onto  the  GaAs  substrate  during  crystal  growth 


kl  =  tvGa]  /  [As4  (g)  ]  ^ 

k2  =  [ASj]  /  (As4(g)]1/4 

k3  =  [AsAs(s>]  [AsGa]  /  [As4(g)]1//4 


(6.7) 


(6.8) 


(6.9) 


From  Eqs.(6.7)  -  (6.9),  yields  the  concentration  of  VGa'  As j ,  and  As^. 


tvGa]  =  kl  tAs4 (g) ] 1'4  /  [Asas(s)J 


-  ^((r-D/4]1/4  =  kx'  (r-1)1/4 
[ASj]  =  k2'(r-l)1/4 
[As^]  -  k3,(r-l)1/2 


(b)  VPE  grown  GaAs  epitaxial  material: 


(6.10) 


(6.11) 


(6.12) 


Arsine  (AsH3)  and  gallium  chloride  (GaCl)  were  used  as  sources  for  As 
and  Ga  in  the  VPE  grown  GaAs.  Let  the  mole  fraction  ratio  of  arsine  and 
gallium  chloride  be  r,  then: 


[AsH3]  /  [GaCl]  =  r 


(6.13) 


If  we  assume  [GaCl]  equal  to  1,  then  [AsH 3]  will  be  equal  to  r.  The 
reaction  of  As4(g)  and  GaCl  can  be  expressed  by 

GaCl  +  (r/4)As4(g)  +  (l/2)H2(g) 


=  GaAs(s)  +  [ (r-1) 1/4]  As4(g)  +  HCl(g) 


(6,14) 


From  Eq.(6.14),  it  is  noted  that  there  is  (r-1)^4  extra  As4(g)  in  the 
reaction  tube,  and  thus  As  atoms  will  be  deposited  on  As  sites  or  As 
interstitial  according  to  Eqs.(6.4)  to  (6.6).  Therefore,  the  concentrations 
of  V^,  ASj ,  and  As^  can  be  expressed  by: 


tVQj]  -  k4[As4 (g)]1/4  /  [ASas(s)]  =  k4'(r-l)1/4 
[ASj]  =  k5' (r-1)1/4 


(6.15) 


(6.16) 

[As^]  *  kg'(r-l)1/2  (6.17) 

From  Eqs.(6.10)  to  (6.12)  and  (6.15)  to  (6.17) ,  it  is  noted  that  for  both 
MOCVD  and  VPE  grown  GaAs,  the  density  of  VGa  is  proportional  to  (r- 
l)1/4;  the  density  of  Asj  is  proportional  to  (r-1)1^4;  and  the  density  of 
As^  is  proportional  to  (r-1)1^. 

6.2.2.  Cooling  process: 

Formation  of  arsenic  anti  site  (As^)  and  arsenic  anti  site-plus-arsenic 
vacancy  (As^V^g)  defects  occurs  when  VGa  migrates  into  As^  or  As^  sites 
under  thermal  nonequilibrium  condition.  Thus,  the  EL2a  level  may  be 
attributed  to  the  As^  defect,  while  EL2b  level  is  due  to  As^V^  .  Fig.6.2 

shows  the  cooling  process  of  GaAs  layers.  This  may  be  further  exemplified 
as  follows.  The  As^  defect  can  be  formed  by  vGa  and  Asj,  which  is  written 
as: 

ASi  +  VQa  =  A (6.18) 

The  ^Ga^s  can  be  formed  by  vGa  migration  to  its  neighbor  (AsAs)  site,  and 
is  given  by: 

ASas<3>  +  =  As^V^  (6.19) 

The  mass  action  law  for  the  reaction  equations  given  in  Eqs.(6.18)  and 
(6.19)  is  described  by: 

Kg  *  [ASgg]  /  [ASj ]  [VQg]  (6.20) 


*b  9  ^GaW  /  UVqJ  [AsAs(s))} 


(6.21) 
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The  ASgg  and  As^V^  defect  concentrations  are  given  respectively  by 

[Mgal  ■  K-,  [ASj]  [V^,]  =  ^"(r-l)1/2  (6.22) 

1*WAS1  -  Kb  [Vg,)  >  Kb,,(c-l)1'/4  (6.23) 

Fig.6.3  shows  the  densities  of  EL2a  and  EL2b  trap  as  a  function  of  [As]/ 
[Ga]  mole  fraction  ratio  as  predicted  by  Eqs.(6.22),  (6.23)  and  the 

published  experimental  data  [75,87,101]^  The  resuits  show  that  for  r  >  1, 
the  density  of  EL2a  trap  is  proportional  to  (r-1)*/2,  and  the  density  of 
EL2b  trap  varies  with  (r-1)^4,  as  predicted  by  Eqs.(6.22)  and  (6.23). 

6.2.3.  Annealing  process: 

The  EL2a  level  may  be  formed  from  EL2b  via  thermal  annealing.  Fig.6.4 
shows  the  formation  mechanisms  of  EL2b  and  EL 2a.  Note  that  EL2b  level  may 
gain  sufficient  thermal  energy,  and  decomposes  into  VGa  according  to  the 
inverse  reaction  of  Eq.(6.19)  when  the  annealing  temperature  is  higher  than 
500  °C  ASgg  is  formed  when  the  migrate  to  the  ASj  site,  according  to 
Eq.(6.18).  This  is  supported  by  the  experimental  data  of  Taniguchi  et 
al^4^  and  Day  et  al^4^. 

6.3.  Detenu  nation  of  Potential  Well  far  the  EL2  Electron  Trap 
from  Field  Enhanced  Rnission  Rate  Analysis. 

As  discussed  above,  the  most  probable  physical  origin  for  the  EL2a 
trap  is  due  to  arsenic  antisite  defect,  and  the  possible  origin  of  EL2b 
trap  is  attributed  to  arsenic  anti  site-plus-arsenic  vacancy  complex  defect. 
The  charge  state  for  these  electron  traps  may  be  revealed  if  the  type  of 
potential  well  is  known.  This  can  be  done  by  studying  the  field  dependent 
emission  rates  deduced  from  the  nonexponential  DLTS  response  for  these 
electron  traps,  as  was  described  in  chapter  IV. 


Fig. 6. 3  The  calculated  density  of  EL2-a  and  EL2-b  electron  traps 
in  GaAs  vs  [As]/[Ga]  mole  fraction  ratios,  as  compared 
with  the  published  data.  For  EL2-a,  Nt  varies  with 


Fig. 6. 4  The  mechanism  of  annihalating  an  EL2-b  trap  and  the  creation 

an  EL2-a  electron  trap. in  GaAs  during  thermal  annealing  process 


The  type  of  potential  well  for  the  EL2a  electron  trap  was  analyzed 
from  the  theoretical  calculations  of  nonexponential  DLTS  response  in  GaAs. 
A  comparison  of  the  calculated  DLTS  response  with  the  measured  DLTS  data 
will  enable  us  to  determine  the  potential  well  for  this  electron  trap.  We 
have  calculated  the  nonexponential  DLTS  response  for  four  different  types 
of  potential  wells,  namely,  the  Coulombic  well,  Dirac  well,  square  well  and 
polarization  well  for  the  EL2a  trap,  and  the  results  are  shown  in  Fig.6.5. 
Fig. 6.6  shows  the  calculated  DLTS  signal  and  the  experimental  DLTS  data  for 
the  Coulombic  potential  well  with  single  and  double  charge  states.  The 
capture  cross  section  (  6^)  of  8xl0”14  cm2  was  used  in  the  calculations.  A 
comparison  of  the  calculated  nonexponential  DLTS  response  with  the  measured 
DLTS  data  shows  that  the  observed  EL2a  trap  is  best  fitted  to  aCoulombic 
potential  well  with  double  charge  state.  Fig. 6. 7  shows  the  theoretical 
nonexponential  DLTS  response  of  Coulombic  potential  well  with  double  charge 
state  for  different  emission  constants  [1/X  =  34.4,  86.6,  172,  344,  and  866 
s-1] . 

The  field  dependent  emission  rates  for  the  Coulombic  potential  well 
can  be  derived  by  taking  into  account  the  three-dimensional  Poole-Frenkel 
effect  and  phonon-assisted  tunneling  effect,  and  can  be  expressed  by: 

a*  H*  .  jn  wt 

enHc/en0=  t  J  f  sin(0)exp(  AEtiAT)d0  d ft/  /  sin(0)  d«  d$>/47r  } 

+/  ***  exp{z-z3/2  [4  (2m*) 1//2  (kT)  3//2/3qrtF] 

0 

(1-  (  4Eti/zkT)5/3]}  dz  (6.24) 

where  e^  is  the  emission  rate  at  zero  electric  field.  2*Etj  is  the  Poole- 
Frenkel  barrier  due  to  the  external  electric  field.  =  q[qFcos(0)  /jj 

io&a]^2r  and  F  is  the  applied  electric  field.  The  first  term  in  Eq.(6.24) 


Fig.6.5  DLTS  scan  of  EL2-a  electron  trap  inGaAs.  Solid  line 
experiment;  dash  line:  theoretical  calculation, 
assuming  Coulombic  well  with  double  charge  state. 


is  due  to  the  three-dimensional  Poole-Frenkel  effect,  and  the  second  term 
is  due  to  the  phonon-assisted  tunneling  effect.  Since  the  electric  field 
varies  with  position  within  the  depletion  region  of  a  reverse  biased  p-n 
junction,  the  emission  rate  is  not  constant  within  the  same  region.  The 
DLTS  response  due  to  the  nonuni  formed  emission  rates  can  in  general  be 
expressed  by  S(tf  *  Z.  exp(-enit1)  -  2  expHfo  t2). 

A  comparison  of  our  theoretical  calculations  of  the  nonexponential 
DLTS  response  due  to  the  field  dependent  emission  rates  with  the  measured 
DLTS  data  reveals  that  the  potential  well  for  the  EL 2a  electron  trap  is  due 
to  Coulombic  well  with  double  charge  state,  and  hence  the  most  likely 
physical  origin  for  this  electron  trap  is  attributed  to  AsGa++.  This 
conclusion  is  consistent  with  the  EL2  modelling  based  on  the  kinetic 
equations  and  defect  structure  as  discussed  in  the  previous  section. 

6.4.  Sununry  and  Conclusions 

In  this  chapter,  we  presented  a  new  defect  model  for  the  EL2  electron 
trap  in  the  undoped  GaAs.  EL2  level  is  attributed  to  two  types  of  native 
point  defects.  One  is  designated  as  the  EL2b  (Ec-0.76eV)  electron  trap, 
and  the  other  is  denoted  as  the  EL2a  (Ec-0.83eV)  electron  trap.  The 
physical  origin  of  EL2a  is  attributed  to  arsenic  anti  site  (AsGa)  defect, 
whereas,  the  physical  origin  of  EL2b  level  is  attributed  to  the  arsenic 
anti  site-plus-arsenic  vacancy  (As^V^g) .  Based  on  this  model,  relationship 
between  the  density  of  EL2a  and  EL2b  trap  levels  vs  [As]/[Ga]  ratio  in  the 
MOCVD  and  VPE  grown  GaAs  was  established.  The  result  shows  that  density  of 
EL2a  trap  is  directly  proportional  to  the  mole  fraction  ratio  of  (r-l)1//2, 
but  the  density  of  EL2b  trap  is  proportional  to  the  mole  fraction  ratio  of 
(r-1)^.  This  prediction  is  supported  by  experimental  data  for  the  MOCVD 
and  VPE  grown  GaAs. 
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VII.  STCDY  OF  GROM- IN  DEEP  LEVEL  D&BCfS  VS  GROWTH  PARAMETERS 
IN  VPE,  LEG,  LPE  and  MOCVD  Grown  GaAs. 

Studies  of  grown- in  deep  level  defects  in  GaAs  grown  by  VPE,  LEC,  LPE, 
and  MOCVD  techniques  under  various  growth  conditions  will  be  described  in 
this  chapter. 

7.1.  Study  of  Grown- in  Deep  Level  Defects  vs  Growth  Parameters 

in  VPE  GaAs  Layers. 

Detailed  study  of  the  grown-in  deep  level  defects  vs  substrate 
orientation  and  gas  phase  stoichiometry  in  the  VPE  GaAs  layers  grown  by  a 
novel  Ga/AsClj/l^  reactor  [which  is  under  high  arsenic  pressure  and  Ga-rich 
conditions]  has  been  made,  using  DLTS  and  C-V  methods.  Density  of  electron 
traps  vs  [Ga]/[As]  ratios  of  2/1,  3/1,  4/1,  5/1,  and  6/1  was  determined 
for  the  VPE  GaAs  layers  grown  on  <100>,  <211A>,  and  <211B>  oriented  semi- 
insulating  (S.I.)  Cr-doped  GaAs  substrates.  Two  electron  traps  with 
energies  of  Ec-0.71  eV  (i.e.,  EB4)  and  Ec-0.83  eV  (i.e.,  EL2a)  were 
observed  in  the  samples  studied.  Results  showed  that  density  of  both 
electron  traps  depends  strongly  on  the  [Ga]/[As]  ratio  for  the  <21lA>  and 
<211B>  oriented  samples,  and  less  strongly  for  the  <100>  oriented  samples. 
For  example,  in  the  <211A>  (i.e.,  Ga-rich  face)  oriented  samples;  the 
densi  ty  of  EB4  and  EL2a  electron  traps  was  found  to  decrease  wi  th 
increasing  [Ga]/[As]  ratio,  while  reverse  trend  was  found  for  the  <211B> 
(i.e.,  As-rich  face)  oriented  samples.  As  for  the  <100>  oriented  samples, 
the  density  of  EL2a  level  may  either  increase  or  decrease  with  increasing 
[Ga]/[As]  ratio.  The  physical  origins  of  this  trap  is  attributed  the 
arsenic  anti  site  (ASQa++),  as  was  discussed  in  chapter  VI.  As  for  the  EB4 
level,  the  results  are  less  conclusive,  and  further  study  is  needed. 
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7.1.1.  Introduction: 


although  the  AsCl^  epitaxial  growth  system  is  the  most  widely  used 
method  for  the  growth  of  microwave  device  materials,  only  limited 
information  is  available  on  how  variation  of  the  [Ga]/[As]  ratio  in  the  gas 
phase  composition  can  affect  the  materials  produced.  In  most  AsCl-j 
reactors,  this  ratio  is  not  easily  controlled  being  a  function  of  source 
efficiency,  wall  deposit  extensiveness,  flow  rates,  etc.  If  stoichiometry 
related  defects  or  defect  complexes  are  important  in  the  VPE  GaAs  growth, 
it  might  be  expected  that  the  [Ga]/[AsCl-j]  ratio  in  the  growth  environment 
is  important.  Miller  et  al  ^86]  have  studied  the  effect  of  gas-phase 
stoichiometry  on  deep  levels  in  the  VPE  GaAs  grown  by  the  hydride  system 
with  [AsH3]/[GaCl]  ratio  varying  from  1/3  to  3/1.  A  deep- level  electron 
trap  (EL2a)  with  energy  of  Ec-0.83  eV  was  observed  in  their  study;  its 
density  was  found  to  depend  on  the  mole  fraction  ratio  of  Ga  to  As  atom. 
Similar  studies  were  also  reported  by  Merenda  ^102J  in  the  VPE  GaAs  grown 
by  the  chloride  system  and  the  MOCVD  grown  GaAs  None  of  these 
studies  have  looked  into  the  effects  of  varying  substrate  orientation  and 
[Ga]/[As]  ratio  on  the  deep- level  defects  in  the  VPE  GaAs  epi layers  as  will 
be  reported  in  this  chapter.  To  facilitate  our  study,  a  novel  AsCl^  vapor 
phase  epitaxial  reactor  [31,103]  was  use<j  t0  grow  the  GaAs  epitaxial 
layers;  the  system  is  capable  of  operation  with  a  variable  controlled 
[Ga]/[As]  ratio.  The  GaAs  epi layers  were  grown  on  three  different 
orientations  (i.e.,  <100>,  <21lA>,  and  <211B»  substrates.  DLTS  and  C-V 
methods  were  employed  to  determine  the  density  of  deep-level  defects  vs 
[Ga]/[As]  ratio  in  the  VPE  GaAs  layer. 
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7.1.2.  Experimental  Details 

j  The  VPE  GaAs  epitaxial  layers  used  in  this  study  were  grown  by  a  novel 
AsClj  reactor  system  131,103]^  The  basic  design  of  this  reactor  is  like  the 
usual  two  bubbler  AsClj  reactor  used  for  the  MESFET  growth  The  novel 

feature  of  our  reactor  is  the  addition  of  a  third  bubbler  which  feeds  a 
cracking  furnace  and  then  a  second  source  chamber  to  allow  operation  in  a 
Gael  rich  mode.  It  uses  a  six-zone  clamshell  furnace  which  rolls  off 
sideways  allowing  an  end  to  end  gas  flow  system.  The  third  or  control 
bubbler  feeds  a  cracking  furnace  that  is  maintained  at  900°C  when  in  use. 

The  GaAs  epi layers  were  grown  on  <100>,  <211A>,  and  <211B>  oriented 
S.I.  Cr-doped  GaAs  substrates  with  growth  temperatures  ranging  from  710°C 
to  750°C.  To  facilitate  our  study  of  the  effect  of  orientation  on  the 
grown-in  defects  in  the  VPE  GaAsf  seven  growth  runs  were  made  for  this 
study  (e.g.,  77A,  77B,  and  77BB  were  grown  in  the  same  run,  etc.).  The  gas 
phase  stoichiometry  was  controlled  by  varying  the  [Ga]/[As]  ratios  from  2/1 
to  6/1.  Al-GaAs  Schottky  barrier  structures  were  fabricated  from  these  VPE 
GaAs  epi  layers  for  the  DLTS  and  C-V  measurements.  Density  and  energy 
levels  of  the  grown-in  defects  in  these  VPE  GaAs  samples  were  determined 
from  the  C-V  and  DLTS  data.  The  results  are  presented  next. 

7.1.3.  Results  and  Discussions. 

In  this  section,  we  present  the  results  of  our  DLTS  and  c-V 
measurements  on  the  VPE  GaAs  samples  described  in  preceding  section.  The 
DLTS  scans  of  electron  traps  as  a  function  of  [Ga]/[As]  ratio  for  different 
runs  and  substrate  orientations  are  shown  in  Fig. 7.1  through  Fig. 7. 9. 
Fig.7.1  shows  the  DLTS  scans  of  electron  traps  as  a  function  of  the 
{Ga]/[As]  ratio,  for  samples  with  <211A>  orientation.  In  general,  one  to 
two  electron  traps  (i.e.,  EB4  =  Ec-0.71  eV,  and  EL2a  =  Ec-0.83  eV)  were 


DLTS  scans  of  electron  traps  for  V 
oriented  substrates  with  [Ga]/[As] 


observed  in  these  samples.  Note  that  sample  37HB  with  the  highest 
[Ga]/[As]  ratio  (i.e.,  6/1)  has  the  lowest  density  for  EL2a  electron  trap, 
while  sample  35  CBwith  the  lowest  [Ga]/[As]  ratio  (i.e.,  2/1)  has  the 
highest  density  of  EL2  electron  trap.  Thus,  for  the  <211A>  orientated 
samples  (i.e„  Ga-rich  face),  the  density  for  EL2a  level  decreases  with 
increasing  [Ga]/(As]  ratio;  this  same  trend  was  also  observed  in  EB4 
electron  trap.  Fig.7.2  shows  the  DLTS  scans  of  electron  traps  as  a 
function  of  the  [Ga]/[As]  ratio  for  epi layers  grown  on  the  <100>  oriented 
substrates.  Results  showed  that  the  density  of  EL2a  level  may  either 
increase  or  decrease  with  increasing  [Ga]/[As]  ratio,  depending  on  the 
concentrations  of  VGa  or  ASj  formed  during  crystal  growth.  To  explain 
this,  modelling  presented  in  chapter  VI  appears  to  be  adequate  for  this 
high  arsenic  pressure  and  Ga-rich  condition.  We  suggested  that  EL2a  level 
is  attributed  to  the  As^  arsenic-anti  site  defect  which  is  formed  by  the 
migration  of  vGa  defect  to  Asj  site  after  crystal  cooling,  as  shown  in 
Eq. (6.18).  VPE  GaAs  layers  are  grown  under  high  arsenic  pressure 
condition.  In  the  case  of  <211A>  oriented  samples,  VGa  should  be 
suppressed  by  the  following  reaction  under  the  Ga-rich  condition. 

Vca  +  Gaj  =  Ga^  (7.1) 

Thus,  the  density  of  EL2a  level  will  decrease  with  increasing  [Ga]/[As] 
ratio,  as  is  the  case  of  <211A>  samples  shown  in  Fig. 7.1.  On  the  other 
hand,  in  the  case  of  <21lB>  oriented  samples,  arsenic  face  will  prevent  the 
reaction  of  Eq.(7J.)  in  the  Ga-rich  case.  At  this  case,  the  concentrations 
of  ASj  and  vGa  levels  shall  be  enhanced  with  increasing  [Ga]/[As]  ratio. 
Therefore,  the  density  of  EL2a  level  will  increase  with  increasing 
[Ga]/[As]  ratio;  this  is  the  case  for  <211B>  oriented  samples.  This  defect 
model  seems  adequate  for  predicting  the  orientation  dependence  of  the 
density  of  EL2a  vs  [Ga]/{As]  ratio,  as  is  shown  in  Fig.  7.9.  Fig.  7.3 
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Fig. 7. 3  DLTS  scans  of  electron  trap  for  VPE  GaA 
(100)  oriented  substrates  with  [Ga]/[As] 


shows  the  DLTS  scans  of  electron  traps  for  sample  35CA  <100>  and  sample 
35CB  <211A>  with  [Ga]/[As]  ratio  equal  to  2/1.  These  two  samples  were 
grown  in  the  same  run  side  by  side  at  715°C.  The  DLTS  results  showed  that 
density  of  EL2a  electron . trap  for  the  <211A>  sample  is  significantly  higher 
than  that  of  the  <100>  sample  (i.e.,  1.24x10**  cm~^  vs.  7.6  x  10*^  cm-"*). 
The  reason  for  this  may  be  attributed  to  the  fact  that  during  crystal 
growth  in  high  arsenic  pressure  condition,  the  density  of  VGa  is  higher  in 
the  <211A>  face  than  that  in  the  <100>  face  for  [Ga]/[As]  ratio  equal  to 
2/1.  The  EB4  electron  trap  was  not  observed  in  both  samples  shown  in 
Fig.7.3.  Fig.7.4  shows  the  DLTS  scans  of  electron  trap  for  sample  75A 
<100>  and  sample  75B  <211A>  with  [Ga]/[As]  ratio  equal  to  3/1.  Note  that 
the  EL2a  level  was  observed  only  in  sample  75A  while  EB4  level  was  observed 
only  in  75B  sample.  The  growth  temperature  for  both  samples  was  held  at 
725°C.  Fig.7.5  shows  the  DLTS  scans  of  electron  traps  for  samples  76A 
<100>  and  76B  <211A>.  The  results  showed  that  both  EL2a  and  EB4  electron 
traps  were  observed  in  these  two  samples.  Density  of  EL 2a  and  EB4  levels 
was  higher  for  sample  76B  <211A>  than  sample  76A  <100>.  Fig.7.6  shows  the 
DLTS  scans  of  electron  traps  for  samples  77B  <211A>,  77A  <100>,  and  77BB 
<211B>.  These  three  samples  were  grown  side  by  side  at  750°C  and  with 
[Ga]/[As]  ratio  equal  to  4/1.  Both  EL2a  and  EB4  electron  traps  were 
observed  in  these  three  samples;  the  density  of  EL2a  and  EB4  levels  was 
highest  for  sample  77BB,  followed  by  samples  77B,  and  77A.  Fig.7.7  show 
the  DLTS  scans  of  electron  traps  for  samples  78B  <211A>  and  78BB  <211B> 
grown  at  725°C  and  with  {Gaj/[As]  ratio  equal  to  5/1.  The  results  showed 
that  EL2a  level  was  the  dominant  electron  trap  for  sample  77B8  while  only 
EB4  level  was  observed  in  sample  78B.  Fig. 7.8  shows  the  DLTS  scans  of 
electron  traps  for  samples  37HA  <100>  and  37HB  <21lA>  grown  at  725°C  and 
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Fig. 7. 5  DLTSscansof electron  trapsinVPE  GaAs  grown  on  (100) 
(211  A)  substrates  with  [Ga]/[As]  ratio  equal  to  4/1. 
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Ga/As  =  4/1 
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electron  traps  for  VPE  Ga As  grown  on  (211  A) ,  (211B)  and 
»nted  substrates  with  [Ga]/[As]  ratio  equal  to  4/1. 
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with  [Ga]/[As]  ratio  equal  to  6/1.  The  results  showed  that  both  EL2a  and 
ES4  level  was  observed  in  sample  37HB.  To  sum  up  the  results  shown  in  Fig. 
7.1  through  Fig.7.8  concerning  the  effects  of  orientation,  [Ga]/[As]  ratio, 
and  the  growth  temperature  on  the  density  of  EL 2a  and  EB4  levels,  we  plot 
in  Fig.7.9  the  density  of  EL2a  level  vs  [Ga]/[As]  ratio  for  three  substrate 
orientations.  In  Fig.  7.9,  it  is  clearly  shown  that  density  of  EL2a  level 
was  found  to  decrease  with  increasing  [Ga)/[As]  ratio  for  the  <211A> 
oriented  samples  and  to  increase  with  increasing  [Ga]/[As]  ratio  for  the 
<211B>  oriented  samples;  for  the  <100>  oriented  samples  the  density  of  EL2a 
level  may  increase  or  decrease  with  increasing  [Ga]/[As]  ratio,  depending 
on  the  vGa  or  Asj  defects  formed  during  crystal  growth.  Table  7.1 
summarizes  the  defect  parameters  deduced  from  our  DLTS  and  C-V  measurements 
on  twelve  VPE  GaAs  samples  prepared  under  different  growth  conditions  and 
substrate  orientations. 


7.1.4  Swzy  and  Conclusions. 

A  study  of  the  grown- in  deep  level  defects  in  VPE  GaAs  vs.  substrate 
orientations,  [Ga]/[As]  ratios,  and  growth  temperature  has  been  carried  out 
in  this  chapter.  One  to  two  electron  traps  (i.e.,  EB4  =  Ec-0.71  eV,  and 
EL2  =  Ec-0.83  eV)  have  been  observed  in  these  VPE  GaAs  samples.  From  the 
results  of  our  DLTS  and  C-V  analysis,  it  is  found  that: 

(1)  background  carrier  concentrations  are  two  to  five  times  higher 
for  the  <211A>  samples  than  the  <100>  samples; 

(2)  defect  density  was  in  general  higher  for  the  <211>  samples  than 
the  <100>  oriented-samples; 

(3)  density  of  EL2a  level  was  found  to  decrease  with  increasing  Ga]/[As] 
ratio  for  the  <211A>  samples,  and  increases  with  [Ga]/[As]  ratio  for 
the  <211B>  samples; 
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(4)  for  the  <100>  samples,  the  density  of  EL2a  level  may  increase  or 
decrease  with  [Ga]/[As]  ratio  depending  on  the  As^  or  VGa 
defect  density  formed  during  crystal  growth; 

(5)  the  EL2a  level  is  attributed  to  the  arsenic  antisite  defect  (AsGa++), 
formed  after  crystal  cooling;  (6)  for  the  EB4  electron  trap,  which  is 
not  commonly  observed  in  the  VPE  GaAs,  and  its  origin  is  not  known. 


Table  7.1  Electron  trap  parameters  in  VPE  GaAs  grown  with  different 
substrate  orientations  and  growth  temperatures. 
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7.2.  Study  of  Grown- in  Deep  Level  Defects 
vs  Growth  Parameters  in  LBC  Grown  GaAs 

In  this  section,  we  present  the  results  of  our  study  of  the  native 
defects  observed  in  the  LBC  grown  bulk  GaAs.  Both  p-  type  GaAs  and  n-  type 
GaAs  have  been  studied.  An  acceptor  level  located  at  80  meV  from  the 
valence  band  edge  has  been  observed  in  the  LBC  grown  p-GaAs.  For  n-GaAs,  a 
study  of  the  effect  of  heat  treatment  on  the  native  defects  on  samples 
annealed  at  200,  300,  and  500°C  for  one  hour  in  hydrogen  (H2)  ambient  has 
also  been  carried  out  in  this  work. 

7.2.1  Intrinsic  Double  Acceptor  Level  in  LBC  Grown  p-GaAs 

Two  p-type  GaAs  specimens  grown  by  the  LBC  method  were  chosen  for  this 
study.  Ohmic  contacts  were  made  by  evaporation  of  Ag-Mn  alloy,  and  an  Al- 
Schottky  barrier  structure  was  used  for  C-V  and  DLTS  measurements. 

Fig. 7. 10  shows  the  DLTS  scan  of  hole  trap  observed  in  both  sample  A 
and  B.  The  apparent  activation  energy  was  obtained  from  the  Arrehenius 
plot  of  the  hole  emission  rate  e^/T  vs  1/kT.  The  measured  activation 
energy  for  this  hole  trap  is  130  meV  from  the  valence  band.  The  true 
activation  energy  of  this  hole  trap  is  80  meV  if  correction  due  to  the 
temperature  dependent  capture  cross  section  is  tak  n  into  account  (i.e., 
tf"p  =  7.1xl0“15  exp(-  AE^/kT)  cm^,  where  AEh  =  50  meV).  Thus,  it  is 
evident  that  this  hole  trap  observed  by  the  DLTS  experiment  is  the  same 
center  observed  in  the  1.441  eV  emission  under  photo  1  urn i nescence 
experimental  iexcitation  and  Hall  measurements.^7^  This  hole  trap  is  a 
double  charged  acceptor  center,  and  is  due  to  the  gallium  anti  site  defect, 

_  _  1C  1  C 

Ga Ag  .  The  density  of  this  hole  trap  was  found  to  vary  from  10XJ  to  10xo 


attributed  to  the  Gaa<=  anti  site  defect. 


-rmally  stimulated  capacitance  (TSCAP)  scans  for  the 
■grown  and  thermally  annealed  LEC  grown  p-GaAs  specimen 


7.2.2  Study  of  Deep  Level  Defects  vs  Annealing  Tesperature 
in  LBC  Grown  n-GaAs 

In  the  following,  the  results  of  the  TSCAP  and  DLTS  measurements  on 
the  LEC  n-GaAs  materials  annealed  in  hydrogen  ambient  are  discussed.  Fig. 
7.11  shows  the  TSCAP  scans  for  the  Schottky  diodes  fabricated  from  LEC 
grown  n-type  GaAs  materials  at  a  reverse  bias  of  -2  volts.  From  this 
figure,  it  is  noted  that  the  capacitance  values  for  samples  annealed  at  200 
and  300°C  for  one  hour  in  H2  ambient  show  little  change  from  the  unannealed 
sample.  However,  as  the  annealing  temperature  raises  to  500°C,  the 
capacitance  value  decreases  by  more  than  50%  from  the  unannealed  sample, 
indicating  a  large  reduction  in  free  carrier  density  in  the  500°C  annealed 
sample.  This  may  be  due  to  the  introduction  of  the  EL2a  electron  trap  by 
the  500°C  annealing  process,  as  is  evidenced  by  our  DLTS  data.  The 
capacitance  steps  observed  in  the  TSCAP  scans  such  as  those  occurred  at  130 
and  240°K  for  sample  annealed  at  500°C  are  due  to  the  electron  emission 
from  the  corresponding  electron  traps,  the  exact  energy  level  for  these 
electron  traps  can  be  determined  from  the  DLTS  technique. 

The  DLTS  scans  of  electron  traps  for  the  LEC  grown  GaAs  annealed  at 
500°C  in  H2  ambient  for  one  hour  are  shown  in  Fig.7.12  and  7.13.  The 
defect  density  deduced  from  the  DLTS  measurements  are  summarized  in  table 
7.2.  For  samples  unannealed  or  annealed  at  200  and  300°C  for  one  hour, 
there  are  three  electron  traps  with  energies  of  Ec-0.35eV,  Ec-0.61eV,  and 
EL2b  *  Ec-0.76eV  oberved  in  these  samples.  As  the  annealing  temperature 
increases  to  500°C,  a  new  electron  trap  with  energy  of  Ec-0.83eV  (EL2a 
level)  is  emerged  into  the  DLTS  scans  along  with  two  other  electron  traps 
(ie.,  Ec-0.35  and  0.61eV),  while  EL2b  level  disappeared;  this  is  shown  in 
Fig.7.13.  The  Ec-0.35eV  level  is  the  main  electron  trap  observed  in  these 


Fig. 7. 12  DLTS  scans  of  electron  traps  in  the  LEC  grown  n-GaAs  specimen. 
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Annealing  performed  in  K 2  ambient  for  one  hour. 
LZC :  Liquid  Encapsulated  Czochralski  technique. 
N,j  measured  at  the  2  um  depletion  width. 


LEC  grown  GaAs  samples.  This  trap  has  also  been  reported  by  Martin  et 
al^77^  in  the  bulk  GaAs  and  by  Fairman  et  al^0^  in  the  VPE  grown  GaAs. 
The  reduction  in  free  carrier  density  with  increasing  annealing  temperature 
along  with  large  capture  cross  section  for  these  electron  traps  indicating 
that  they  are  acceptor  type  electron  traps.  Another  interesting  study  on 
these  traps  is  concerned  with  the  spatial  dependence  of  the  defect  density 
for  each  electron  trap.  This  is  obtained  by  performing  the  bias  dependence 
of  the  DLTS  scans,  and  the  results  are  illustrated  in  Fig.7.14.  From  the 
DLTS  scans  as  a  function  of  the  reverse  bias  voltage,  the  spatial 
dependence  of  trap  density  is  determined  for  each  trap.  Fig.  7.15  show  the 
spatial  dependence  of  the  electron  traps  observed  in  the  200°C  annealed 
samples.  It  is  shown  that  for  all  the  samples  studied  here,  the  defect 
density  is  highest  at  the  surface,  and  decreases  as  it  moves  away  from  the 
interface  of  the  metal  GaAs  Schottky  contact;  the  trap  density  becomes 
constant  deep  into  the  bulk  region  (ie.,  2um)  of  the  GaAs  substrate, 
indicating  that  all  these  electron  traps  are  bulk  related  defects.  Mircea 
et  al^®l  have  also  observed  a  similar  defect  density  profile  for  the  Ec- 
0.61eV  level  in  which  the  density  of  Ec-0.61eV  was  increased  by  the 
mechanical  damage  near  the  surface  of  GaAs  specimen. 

Fig.  7.16  shows  the  spatial  dependence  of  the  Ec-0.35eV  electron  trap 
in  the  unannealed  and  annealed  samples.  The  results  showed  that  density  of 
the  Ec-0.35eV  trap  was  decreased  with  increasing  annealing  temperature. 
EL2b  trap  was  annealed  out  followed  by  a  500°C  annealing  in  H2  ambient  for 
one  hour.  Recently,  Day  et  al^4]  has  observed  the  Ec-0.35eV  and  EL2b  = 
Ec-0.76eV  traps  in  the  MBE  grown  n-GaAs,  and  found  that  both  Ec-0.35eV  and 
EL2b  traps  can  be  annealed  out  at  800°C  for  1/2  hour.  To  explain  their 
results,  the  EL2  model  proposed  in  previous  chapter  appears  to  be  adequate 
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Fig. 7. 15  Density  of  electron  trapsinthe  200  C  annealed  GaAs 
specimen  as  a  function  of  the  depletion  layer  width. 
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y  ec~0»35  eV  electron  trap  vs  depletion  layer  width  in 
n-GaAs  for  different  annealing  temperatures. 


for  this  case.  The  physical  origin  of  EL2b  trap  is  attributed  to  AsGaVfts 
complex.  At  500°C  annealing  temperature,  the  EL2b  level  (i.e.,  the  As^V^g) 
will  decompose  into  vGa'  and  its  reaction  equation  is  given  by: 


AsGaVAs  "  ^As  +  vGa 


(7.2) 


Note  that  at  high  annealing  temperatures  (eg.,  Ta  >  500°C),  the  trap  may  be 
formed  by  the  migration  of  vGa  defect  to  the  ASj  site  to  form  an  arsenio- 
antisite  defect  as  was  predicted  by  Eq.(6.18).  This  defect  model  seems 
adequade  for  explaining  change  of  the  EL2b  trap  to  EL2a  level  at  high 
annealing  temperatures  as  shown  above. 


7.3.  GcowD-ln  Deep  Level  Defects  vs  Growth  Parameters 
in  the  LPE  n-GaAs  Layers. 

Samples  of  LPE  GaAs  were  prepared  for  studying  grown-in  defects 
vs  growth  temperature  and  growth  rate  in  the  high  purity  n-GaAs 
epi layers.  These  samples  were  prepared  by  infinite  solution  melt  LPE 
technique  grown  at  700  and  800°C  with  cooling  rate  of  0.4  and 
l°C/min.,  respectively.  Fig. 7.17  shows  the  DLTS  scans  of  electron 
trap  (EB4=Ec-0.71eV)  for  samples  grown  at  700  and  800°C  with  cooling 
rate  of  l°C/min..  Note  that  for  the  same  cooling  rate,  sample  grown 
at  800°C  has  defect  density  six  times  lower  than  that  of  the  700°C 
grown  sample.  On  the  other  hand,  for  sample  grown  at  700  and  800°C 
with  a  cooling  rate  of  0.4°C/min.,  no  electron  traps  were  detected  in 
these  samples.  The  results  of  our  DLTS  and  C-V  measurements  on  the 
LPE  GaAs  were  summarized  in  table  7.3.  The  background  concentration 
in  these  samples  is  around  10 16  cm-3. 


Fig. 7. 17  DLTS  scans  of  electron  trap  for  two  LPE  GaAs  grown  at 

700  and  800  C  with  1  C/min.  cooling  rate. 


7.4.  Growo-in  Deep  Level  Defects  vs  Growth  Parameters 
in  MOCVD  Grown  n-GaAs  Epi- layers. 

Characterization  of  native  defects  in  n-GaAs  epi layers  grown  on 
S.I.  GaAs  and  S.I.  Ge  substrates  by  the  MOCVD  technique  has  been 
studied.  The  Sn-doped  GaAs  epi  layer  on  S.I.  GaAs  substrate  was  grown 
at  700°C  for  30  min.  with  [Ga]/{As]  ratio  of  11/1,  and  the  dopant 
density  was  4.1x10^  cm-3.  As  for  the  Sn-doped  GaAs  epi  layer  on  the 
S.I.  Ge  substrate,  the  sample  was  prepared  at  706°C  for  30  min.  with 
[Ga]/[As]  ratio  of  7/1;  the  background  concentration  was  found  to  vary 
between  10^®  to  1.9x10^  cm-3.  Pig. 7. 18  shows  the  DLTS  scan  for  n- 
GaAs  layers  on  S.I.  GaAs  substrate.  The  result  shows  that  there  is 
one  electron  trap  with  energy  of  Ec-0.6eV.  This  trap  level  is 
corresponding  to  the  Ec-0.61eV  for  the  LEC  grown  GaAs  samples.  As  for 
the  GaAs  layer  on  S.I.  Ge-substrate,  one  hole  trap  with  energy  of  Ev+0.75eV 
and  the  density  of  1.7xl013  cm-3  was  observed  in  this  sample. 

7.5.  Sumaary  and  Conclusions. 

From  the  results  of  our  study  of  the  grown-in  defects  in  the  LEC,  VPE, 
LPE,  and  MOCVD  GaAs  samples,  it  is  clear  that  grown-in  defects  and  the 
background  concentration  are  sensitive  to  the  growth  conditions, 
stoichiometry,  and  substate  orientations.  Table  7.4  shows  the  trap  level 
observed  in  this  study.  The  main  conclusions  are  listed  as  follows: 

(1)  EL2a  =  Ec-0.83eV  electron  trap  is  attributed  to  the  (AsGa++)  arsenic- 
anti  site  defect  which  was  observed  in  VPE  GaAs,  and  in  500°C  annealed 
(H2)  LEC  grown  GaAs. 

(2)  EL2b  a  Ec-0.76eV  electron  trap  is  attributed  to  the  AsGgVAs 
complex  which  was  observed  in  the  LEC  grown  GaAs.  This  level  was 
annealed  out  by  a  500°C  hydrogen  gas  heat  treatment. 
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(3)  Ev+0.08eV  hole  trap  was  observed  in  the  LEG  grown  p-GaAs  samples. 
The  physical  origin  of  this  hole  trap  is  believed  to  be  due  to 
Gafts  gallium-anti  site  defect  as  confirmed  by  the  Hall  effect  and 
photoluminescence  measurements. 

(4)  Reducing  the  temperature  cooling  rate  during  crystal  growth  is 
beneficial  for  reducing  the  density  of  grown-in  defects  in  LPE 
GaAs  epilayes. 


Fig.7.18  DLTS  scan  of  electron  trap  in  MOCVD  GaAs  grown  on  S.I£aAs  substrate 


VIII.  SQMMOr  AND  CONCLUSIONS 

A  detailed  theoretical  and  experimental  study  of  the  grown- in  deep- 
level  defects  and  radiation-induced  defects  in  GaAs  has  been  carried  out  in 
this  research.  The  main  contributions  from  this  research  include:  (i) 
Development  of  a  new  defect  model  for  predicting  native  defects  and  defect 
density  vs  growth  parameters  in  GaAs,  (ii)  modelling  and  new  interpretation 
of  the  EL2  electron  trap  in  GaAs,  (iii)  determination  of  the  energy  level 
and  trap  density  for  a  large  number  of  deep- level  traps  observed  in  GaAs 
grown  by  LEC,  VPE,  LPE,  and  MOCVD  techniques  under  different  growth  and 
annealing  conditions,  (iv)  determination  of  the  potential  well  for  several 
electron  traps  in  GaAs  from  the  field  enhanced  emission  rates  deduced  from 
the  DLTS  data,  which  enable  us  to  determine  the  charge  state  of  each  trap 
in  GaAs.  The  research  findings  described  in  this  report  has  led  to  the 
following  conclusions: 

(1)  A  theoretical  model  for  interpreting  the  nonexponential 
capacitance  transients  due  to  electric  field  dependent  emission  rate  of 
trapped  charge  has  been  developed.  Based  on  the  nonexponential  capacitance 
transient,  theoretical  calculations  of  DLTS  response  for  the  deep  level 
traps  in  GaAs  were  made  using  five  different  potential  wells;  namely,  the 
Coulombic  well,  Dirac  well,  square  well,  polarization  well,  and  dipole 
well.  A  comparison  of  the  theoretical  calculation  of  the  nonexponential 
DLTS  response  with  the  DLTS  data  for  each  trap  level  enable  us  to  determine 
the  potential  well  for  each  trap  in  GaAs.  This  method  has  been  applied  to 
determine  the  charge  state  and  the  physical  origins  of  the  EL2  electron 
trap  in  GaAs. 
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(2)  A  theoretical  model  for  predicting  native  point  defects  in  GaAs 
has  been  developed.  The  chemical  thermodynamic  principles  are  employed  to 
derive  the  density  of  vacancy  and  interstitial  defects  during  crystal 
growth  and  under  thermal  equli librium  condition.  After  crystal  growth,  the 
thermal  kinetic  equations  are  employed  to  predict  the  anti  site  defect.  The 
results  show  that  :  (i)  high  purity  GaAs  material  can  be  grown  for  the  low 
arsenic  pressure  case  under  optimum  cooling  condition,  (ii)  GaAs  grown 
under  higher  arsenic  pressure  condition  (e.g.,  LBC,  VPE  and  MOCVD  grown 
GaAs  )  will  produce  more  native  point  defects  than  under  lower  arsenic 
pressure  condition  (e.g.,  LPE  grown  GaAs),  (iii)  arsenic  anti  site  defect 
can  only  be  observed  in  GaAs  grown  under  As-rich  or  high  arsenic  pressure 
condition;  this  defect  can  not  be  produced  under  low  arsenic  pressure  and 
Ga-rich  conditions. 

(3)  A  new  defect  model  for  the  EL2  electron  trap  in  GaAs  was 
developed  in  this  study.  The  EL2  electron  trap  is  attributed  to  two 
different  native  point  defects.  One  is  designated  as  the  EL2a  (Ec-0.83eV) 
electron  trap,  and  the  other  is  denoted  as  the  EL2b  (Ec-0.76eV)  electron 
trap.  The  physical  origin  for  the  EL2a  level  is  attributed  to  the  arsenic 
antisite  (AsGa++)  defect,  whereas,  the  physical  origin  for  the  EL2b  level 
is  attributed  to  the  arsenic-  antisite  and  arsenic-  vacancy  complex 
(As^v^g) .  Based  on  this  model,  relationship  between  the  density  of  EL2a 
and  EL2b  trap  levels  vs  [As]/[Ga]  mole  fraction  ratio  in  the  MOCVD  and  VPE 
grown  GaAs  is  established.  The  result  shows  that  density  of  EL2a  trap  is 
proportional  to  the  mole  fraction  ratio  of  (r-1)^^;  while  the  density  of 
EL2b  trap  is  proportional  to  the  mole  fraction  ratio  of  (r-l)1//4,  where  r  * 
[As]/[Ga].  This  prediction  is  in  good  agreement  with  experimental  data 
reported  for  the  MOCVD  and  VPE  grown  GaAs. 
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(4)  From  the  results  of  our  DLTS  study  of  the  grown- in  defects  in  the 
LBC,  VPE,  LPE,  and  MOCVD  GaAs  samples,  it  is  shown  that  grown- in  defects 
and  the  background  concentration  are  sensitive  to  the  growth  conditions, 
stoichiometry,  and  substrate  orientations.  The  main  technical  findings 
from  this  study  include:  (i)  EL2a  electron  trap  was  observed  in  VPE  grown 
GaAs  and  in  LEC  grown  GaAs  annealed  at  500°C  in  H2  ambient;  (ii)  EL2b 
electron  trap  was  observed  in  the  unannealed  LEC  grown  GaAs,  and  its 
density  can  be  reduced  by  low  temperature  thermal  annealing;  (iii)  a 
shallow  hole  trap  with  energy  of  Ey  +  0.08  eV  was  observed  in  the  LBC  grown 
p-GaAs.  The  physical  origin  of  this  hole  trap  is  attributed  to  the  gallium 
antisite  (Ga^g)  defect;  (iv)  using  slow  temperature  cooling  rate  after 
crystal  growth  is  beneficial  for  reducing  the  density  of  grown-in  defects 
in  the  LPE  grown  GaAs.  The  theoretical  modelling  of  native  point  defects 
and  EL2  electron  trap  in  GaAs  described  in  this  work  was  found  in  good 
agreement  with  the  experimental  results  for  the  deep- level  defects  in  GaAs 
grown  by  different  growth  techniques  and  growth  conditions  presented  in 
this  report. 

Studies  of  radiation  induced  deep- level  defects  in  the  LPE  grown  GaAs 
irradiated  by  one-MeV  electron  irradiation  and  low  energy  proton 
irradiation  under  different  irradiation  and  annealing  conditions  have  also 
been  carried  out  in  this  research  program,  and  the  results  have  been 
presented  in  the  1983  AFOSR  Annual  Technical  Report. 
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